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MOLECULAR DIFFUSION RATES IN: SUPERCRITICAL 
WATER VAPOR ESTIMATED FROM VISCOSITY DATA 


MATT WALTON 


Department of Geology, Yale University, New Haven, Connecticut 
ABSTRACT. According to the Stokes-Einstein equation the diffusivity of a particle in 
a gas is a function of the radius of the particle, the absolute temperature, and the viscosity 
of the gas, Recent theoretical and experimental results are used to extrapolate values for 
the viscosity of supercritical water vapor up to 800°C, and 2,500 kg cm™ pressure. The 
function Dr, is evaluated in this range, where D is the diffusivity and r, the radius of a 
molecular particle in solution in supercritical water vapor, lt is shown that molecular com- 
ponents of rocks should have diffusion rates on the order of 10° to 10 cm*sec™ over a 
wide range of geologic pressures and temperatures, As an example these results are applied 
to an evaluation of the diffusivity of silica in supercritical water vapor, and it is shown 
that the mass transport of silica by diffusion is potentially a relatively rapid geologic 
process where supercritical water vapor exists at high pressures in the intergranualr spaces 
in rocks, even where the effective directional porosity of the rock is several orders of mag- 
nitude less than one, Values for viscosity of water, Dr,, solubility in water of SiOz in gm 
cm“, mass transport of SiO., and mean square displacement of SiOz in supercritical water 
vapor are given graphically for various temperatures up to 800°C. and pressures up to 
2,500 kg cm™. 

The properties of water in the range of several hundred to a thousand 
degrees centigrade and up to several thousand atmospheres of pressure are of 
great interest, especially in geology. Thanks largely to the work of Kennedy 
(1950; Holser and Kennedy, 1958, 1959) considerable equilibrium PvT data 
are now available in this range. However, the so-called non-equilibrium or 
transport properties have as yet received little systematic attention from geol- 
ogists. These properties, viscosity, thermal conductivity, diffusivity, ete., are 
the properties which determine the capacity of water to serve as a vehicle for 
the transfer of matter and energy. In this paper an attempt is made to estimate 
the capacity of supercritical water to serve as a medium for the diffusion of 
molecular particles. 

Available experimental data on the non-equilibrium properties of water 
do not extend very far into the range of geologic pressures and temperatures. 
To get some idea of how supercritical water behaves in this range it is neces- 
sary to make empirical or theoretical extrapolations. Unfortunately no common 
gas is more intractable to theoretical extrapolation than water vapor because of 
the extreme polarity of the water molecule. 

Diffusivity of a particle in a fluid medium is given by the well-known 
Stokes-Einstein equation: 

tT 
(1) D 


Where D is diffusivity, R is the gas constant, T the absolute temperature, 
Avogadro's number, r, the radius of the particle, and » the viscosity of the fluid 
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medium. It is assumed that the concentration of the diffusing particle is too low 
to affect the properties of the fluid, that the particle is essentially spherical and 
large enough so that the fluid behaves as a continuous medium, and that there 
is no interaction between the particle and the fluid other than the transfer of 
momentum, Since these assumptions are not true with respect to molecular 
particles in a real gas, the equation is at best an approximation, It has been 
found, for example, in gases at high pressure that better agreement with experi- 
mental results is obtained if the numerical coefficient in the denominator is 4 
instead of 6 (Prof. P. A. Lyons, Chemistry Department, Yale University, per- 
sonal communication), On the other hand, the effect of non-ideality of a real 
gas is to some degree manifested in this equation by variation in » with T and 
P, so that if a good approximation can be made of », the departure from 
ideality may be partly compensated. 

I propose to use this equation to get an approximate notion of the rates 
at which molecular components of rocks and minerals might diffuse through 
supercritical water vapor, If the results are of even order of magnitude validity 
they will be highly significant for geology. 

Evaluating the constant (R/4eN) in equation (1) and transposing: 


(2) Dr, 1.1 10-“erg 
The viscosity is, of course, a specific property which varies with T and P, The 
problem is to extrapolate experimental values for viscosity into the range of 
geologic temperatures and pressures. Recently Fano. Hubbell, and Beckett 
(1956) derived functions for the viscosity of water at higher temperatures and 
pressures which allow in some measure for the non-ideality of water vapor by 
adjusting the Chapman-Enskog theory of dense gases (Hirschfelder, Curtiss, 
and Bird, 1956, p. 649) to the PvT data of Kennedy. Their calculated values 
for viscosity are almost all within 10 percent and generally within 5 percent of 
the measured values of Timroth (1940), Sigwart (1936a.b), and Jackson 
(1949). which extend to 600°C and several hundred atmospheres pressure. 

Fano, et al. (1956, p. 11) give an empirical function for the so-called 
zero pressure viscosity of water vapor (which is for all practical purposes the 
viscosity at | atmosphere within the geologic range of temperature) : 

(3) 0.361 T — 10.2, for T < 800°RK. 

9.37 
3.15 + 0.001158 

The values for », derived from these equations compare closely with values 
obtained by Krieger (1951) in an extrapolation of zero pressure viscosity to 
higher temperatures in which an effort was made to take into account the 
dipole interaction of the polar molecule. 

To obtain the viscosity at higher pressures the following function is given: 

(5) n/m 1 + 0.175 bp + 0.8651 b*p? 
in which yy is m1: atmos.) at stated T, p is the density of water vapor at stated 
T and P, and b is obtained from the following: 
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Water Vapor Estimated from Viscosity Data 


(6) b <x 10-7 


in which M is molecular weight of water. Values for the viscosity of water 
Vapor at various temperatures and pressures obtained from these functions are 


VISCOSITY of WATER VAPOR 
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Fig. 1. Viscosity, ” of supercritical water vapor calculated from empirical functions 
of Fano, et al, (1956). 
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presented in column 4 of table 1, and are shown graphically in figure 1, Values 
for density used in the calculation were obtained from Kennedy's tables (Ken- 
nedy, 1950; Holser and Kennedy, 1958, 1959) and are given in column 3 of 
table 1. 

Better than order-of-magnitude significance for these calculated viscosities 
seems reasonable since the values lie between the values for liquid water in 
equilibrium with water vapor just below critical pressure and temperature, 
which are on the order of 10~° poise, and the values for water vapor at low 
poise. The 
viscosity of a liquid tends to decrease with increasing temperature, while the 


pressures and temperatures. which are somewhat less than 10 


viscosity of a gas increases, It is not surprising that the viscosity of supercriti- 
cal water vapor at high pressures exceeds the viscosity of very hot water just 
helow critical temperature. The lowest experimental value obtained by Timroth 
for the viscosity of liquid water is 6.8 * 10~* poise at 360°C, and 200 kg 
em~* pressure. This is less than the viscosity for supercritical water vapor in 
the upper range of pressures for which calculations have been made. 

Values for the function Dr, for supercritical water vapor obtained from 


equation (2 


figure 2. The characteristics of this function have interesting geologic implica- 


) are given in column 5 of table 1, and are shown graphically in 


tions. Once beyond the critical region of pressure and temperature the func- 
tion becomes slowly varying with respect to both T and P over a fairly wide 
and significant range of geologic pressures and temperatures, To the extent 
therefore that the diffusion of molecular particles is dependent on the viscosity 
and kinetic energy of the water vapor, diffusion rates should vary within less 
than an order of magnitude over a considerable geologic range, Moreover, it 
is probable that the kinds of particles, such as molecues and hydrated ions, 
which diffuse in water vapor have radii (r,) ranging from a little less than one 
to three or four angstroms, ( msequently diffusion coetlicients in super ritical 


cm set under 


water vapor should fall in the range of 107° to 10 
conditions ranging from high temperature hydrothermal processes to moderate 
and high-grade metamorphism and magmatic processes, so long as the concen- 
tration of solutes in the Vapor! does not become large enough to seriously modi- 
fy the properties of the gas. 

lo illustrate the application of this function, diffusion rates are calculated 
for silica in water vapor. The values are given in column b of table 1 and 
shown graphically by referring to the lower scale on the abscissa of figure 2. 
It is assumed that silica occurs in solution in supereritical water vapor as a 
tetrahedrally coordinated Si(OH), molecule (Wasserburg, 1958), It would 
appear reasonable to assume that Si(OH), is structurally analogous to other 
hydroxy-acids of the form A(OH), which have been investigated in which the 
OH ions cluster around each cation and adjacent clusters are linked by 
(OH)-(OH) bonds. These bonds range in leneth from 2.70 A to 2.82 A 


(Wells, 1945). The ionic radius of (OH) > is approximately the same as that 
of O?> (Wyckoff, 1931, p. 192). and so I will take as an estimate of the radius 
of the free particle, Si(OH),, the observed interatomic distance in the tetra- 
hedral ion SiO,*~ as given by Pauling (1945, p. 240), which is 1.60 A, plus 
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Fig. 2. Stokes-Einstein diffusion function for supercritical water vapor. Upper scale 
on abcissa gives values for product of diffusivity, D, and radius of solute particle, rs. 
Lower scale gives values for diffusivity of hydrated silica molecule, Dsicom,. 


one half the average (OH)-(OH) bond, which is 1.38 A, giving a value for 
r, of 2.98 A 

The calculated diffusion rate for silica in supercritical water vapor ranges 
from about cm*sec~! to em*sec~'. Diffusion rates of this order are 
large enough to produce very large transfers of material within geologically 
short periods of time. By comparison the diffusivities of common salts in water 
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at s.t.p. are on the order of 10~° to 10~° cm*sec—*, The diffusion rate of 
glycerol at 0.1 molar concentration in water at 10°C is 4.13 X 10~° cm*sec~* 
(molecular weight of glycerol, 92; Si(OH) 4, 96). 

As a basis for estimating the geologic significance of diffusion in super- 
critical water vapor a system may be considered where in one region the vapor 
is saturated with silica at a given temperature and pressure and in another 
region the silica is virtually quantitatively removed by reaction to form an 
insoluble compound, Such a system would be approximated by the experiment 
described by Bowen and Tuttle (1949) where crucibles containing finely di- 
vided silica and MgO are placed side by side in a hydrothermal bomb, In a 
very short time forsterite crystals form in the crucible containing MgO by 
diffusion of silica through the vapor, but no such crystals form in the silica 
crucible because MgO and forsterite are virtually insoluble in supercritical 
water vapor. The concentration gradient in the vapor in such a system should 
approach the solubility of silica in the vapor. 

| have converted the data given by Morey (1957) on the solubility of 
silica in water vapor to grams per cc, using Kennedy’s tables for the specific 
volumes. These data are plotted in figure 3. Values taken from these curves for 
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Fig. 3. Solubility of silica in supercritical water vapor in gm cm™, 
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the concentration of silica in saturated water vapor are given in column 7 of 
table 1. Morey’s data on solubility covers only part of the range of P and T 
for which I have calculated viscosity and diffusivity, but enough to illustrate 
the possible magnitude of silica transport in supercritical water vapor. 

Fick’s equation for mass transport where diffusion is driven by a concen- 
tration gradient is: 

(AC)t 
l 


where m is the mass of solute diffused through a cross-section area, A, in time 
t, D is the diffusivity of the solute in the solution, and AC is the change in 
concentration over a distance /, Assuming the concentration falls from satura- 
tion to zero between two regions at a given T and P, the mass transport is: 


(7) m 


where: 
(9) DC’ 


D being the diffusivity of the solute at stated T and P, and C’ its solubility in 
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Fig. 4. Calculated mass transport of silica by diffusion in supercritical water vapor. 
Upper scale on abcissa gives values for a concentration gradient, AC, of saturation to 0 
in 1 em, Lower scale gives values for a concentration gradient of saturation to 0 in 1 m. 
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grams per ce at stated T and P. Values for M’ for silica in supercritical water 
vapor are given in column 8 of table 1 and shown graphically in figure 4 ac- 
cording to the upper scale on the abscissa. M’ in e.g.s, units is of course the 
mass transport of silica per second through a square centimeter of cross section 
where the concentration drops from saturation to zero in one centimeter. For 
a system more nearly approaching ordinary geologic magnitudes the mass 
transport per year through a square meter with a gradient from saturation to 
zero concentration in one meter is given in column 9 of table 1 and shown 
graphically according to the lower scale on the abscissa of figure 4, It is seen 
that the mass transport ranges up to several kilograms per year, and it is evi- 
dent that where thousands to hundreds of thousands or millions of years are 
available mass transport in supercritical water vapor can be of very great 
geologic magnitude, even with concentration gradients much smaller than have 
heen assumed here. 

Another value which gives some idea of the magnitude of diffusion is the 
so-called mean square displacement, x. 


(10) x" 2Dt 


x is the distance from an interface in a diffusion medium at which the cone-n- 
tration of a solute will reach one half its concentration at the interface in time t. 
Values for x where t | year are given in column 10 of table 1 and shown 
graphically in figure 5. The data indicate that if supercritical steam comes in 
contact with rock materials and becomes saturated with dissolved silica, it will 
reach half saturation at distances of one to several meters from the rock face 
in a year, assuming no motion of the vapor itself; tens of meters in a hundred 
years, thousands of meters in a million years. Since the solubility of silica in 
steam is not large, little mass transport is involved in establishing the half- 
value for concentration at a given distance. What is indicated here is some idea 
of the time necessary to establish a steady concentration gradient in systems of 
various dimensions by diffusion in supercritical water vapor. 

If the diffusivities of molecular particles in high-pressure supercritical 
water vapor are even crudely approximated by the Stokes-Einstein equation 
using the values for viscosity derived here, then it is evident that diffusion in 
supercritical water vapor is potentially a very rapid and large-scale process, as 
geologic processes go, for any substance that has appreciable solubility in 
steam. The Stokes-Finstein equation gives useful results for gaseous diffusion 
in many organic fluids, and the values obtained here for viscosity, as well as 
the diffusivities derived from them, bear a reasonable relationship to the ex- 
perimentally investigated range of measured viscosities in steam and diffusion 
rates in liquid water. Consequently it seems worthwhile to give serious con- 
sideration to the geologic impli ations of these results. 

Up to now I have not dealt with the application of diffusion in super- 
critical water vapor to any realistic geologic model, There may be situations in 
voleanic or hydrothermal vein systems where an extensive body of supercritical 
water vapor exists for a time in an otherwise free space, Surely, however, such 
situations are relatively rare and of short duration. (As a petrologist, | am in- 


clined to look on ore deposits in general as most uncommon rocks). Probably 
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Fig. 5. Mean square displacement of silica diffusing in supercritical water vapor 
from a saturated interface. 


supercritical water vapor most commonly exists as a porefluid or intergranular 
dispersed phase occupying very small free spaces in rocks, The dispersed phase 
may occupy such small spaces that the properties of the fluid are modified by 
adsorbtive forces on enclosing mineral grain surfaces. Nevertheless it is evident 
in many metamorphic terrains that large amounts of water have evolved from 
or been introduced into extensive bodies of rock. It is also evident that in the 
crystallization of magmas at depth, large amounts of water are evolved and re- 
act with the igneous rock and its surroundings, Korzhinsky (1936) points out 
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that in metasomatic processes where chemical transport is taking place 
throughout a body of rock, the proportion of rock that can be affected by the 
flow of solutions in permeable openings is very small: 

. the rock does not react [to a significant extent] with flowing solutions, it 
reacts with an immobile solution permeating the submicroscopical pores of the 
rock and serving as an intermediary between the whole rock mass and the flowing 
solutions. The introduction or leaching of rock components are affected by diffu- 
sion through this immobile subcapillary medium . . . (p. 58, English summary) 


The experimental results obtained by Duffell (1937) tend to support this con- 
cept of solute transport in pore solutions. 

Some rather limited but highly suggestive experiments by Garrels, Dreyer, 
and Howland (1949) indicate that the rate of movement of a substance through 
a rock by diffusion in a pore fluid is nearly independent of the fabric, so long 
as the rock does not interact with the diffusing substance or the properties of 
the fluid become modified by interaction with the pore walls. One can visualize 
by way of analogy an area traversed by a number of paths, The rate of speed 
at which one may walk is more or less fixed, so that the time used by any in- 
dividual to cross the area depends only on the croookedness of the path taken. 
The total volume of traffic, however, in a given time is also a product of the 
number of available paths. 

7), the term A must be 
modified by a factor, a, less than one. which is the fraction of the average cross 


To return to Fick’s equation for mass transport ( 


section occupied by effective diffusion paths. The term / must be modified by 
a factor, j, greater than one, which is the ratio of the average length of the 
diffusion paths to the straight-line distance through the rock for a given change 
in concentration. The equation becomes: 


(AC)t 


(11) D Aa 
] 
Garrels, et al. (1949) use a slightly different notation but in effect 


a 
define — as a property of any given rock which they call the effective direc- 


tional porosity. It is a factor which reduces the cross section area of a rock 
for a given concentration gradient to the effective cross section of the free dif- 
fusion paths through the pore fluid in the rock, as though these paths were 
consolidated into a single hole through the rock. 

To return to the diffusion of silica in supercritical water vapor, the data 
plotted on figure 4 indicate that hundreds to thousands of grams of silica may 
be transported per year through an area | m* with a concentration gradient 
equal to the solubility of silica at the given P and T across a distance of 1 m, 
if the vapor pressure of water exceeds a few hundred atmospheres and the tem- 
perature a few hundred degrees C. This means that at a surface where silica 
is being remorved from solution a layer of pure quartz a centimeter thick could 
be built up over the entire area of diffusion in a few years at commonly at- 
tained geologic pressures and temperatures. The time involved in many geol- 
ogic processes such as metamorphism appears to be on the order of 10° or 
perhaps 10° years. Consequently the product of effective directional porosity 
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and concentration gradient may be reduced by corresponding orders of magni- 
tude and still allow geologically important transport in the intergranular vapor 
phase. 

Garrels, et al. (1949) determined a few values for effective directional 
porosity from measurements of diffusion in water saturated limestones at tem- 
peratures of 25 C to 65 C and got values ranging from about 0.03 to 0.8 with 
some indication of an increase in effective directional porosity with an increase 
in temperature. Duffell found measurable penetration of diffusing ions in 
aqueous solutions in crystalline limestones with porosities of less than one per- 
cent and little or no measurable permeability. These results can only suggest 
that until actual experimental measurements are made of diffusion in inter- 
granular supercritical water vapor in rocks at temperatures and pressures in 
the range of metamorphic, magmatic and hydrothermal processes it would be 
well not to dismiss this as a transport mechanism corresponding in rate and 
magnitude to other major geologic transport processes. 

In fact the basic diffusion rates in the vapor are so rapid that the real rate 
controls in the geologic realm are likely to be factors that limit effective activity 
gradients, or else phenomena such as adsorbtion on boundary surfaces in very 
small openings which modify the properties of the vapor phase, Complex 
coupled reactions may be involved which control the relationship of the rate at 
which the diffusing particle enters the vapor phase to the rate at which it is 
removed elsewhere in the system. The reaction at either interface may produce 
a product layer and diffusivity in the product layer will influence the rate of 
further reaction. Water in intimate contact with the surfaces of silicates tends 
to undergo ion exchange with the surface, hydrolizing the silicate surface and 
taking metal ions into solution (Stevens, 1934) (Frederickson, 1951). Ulti- 
mately a steady state tends to be reached in the pH of the water, which is 
always higher than 7. In capillary and subcapillary pores the base ions in solu- 
tion must tend to become balanced against hydroxyl ions and more or less 
neutral with respect to the molecular forces in the adsorbed water skin on the 
surfaces. This may more than double the diffusion rates of base ions over 
molecular diffusion rates in neutral solutions because the mobility of hydroxyl 
is four or five times at great as the mobility of the ions of the common bases 
in rocks, and hydroxyl “drags” the metal ions along. No doubt there are other 
ways in which activity gradients are affected by factors other than concentra- 
tion. For example, Wang (1954) has discussed an effect which reduces the 
activation energy of diffusion in liquid water owing to the influence on short 
range order in water of small surface charges on large ions or molecules in 


solution, This has an exponential effect on increasing diffusion rates, and may 


he especially important in subcapillary openings where the influence of surface 
charges might be large. Only empirical data can evaluate these effects. 

In the meantime the calculated values given here should help to give some 
quantitative significance to speculations concerning the role of diffusion in 
high-pressure water vapor in geologic processes and perhaps serve as a guide 
to the design of appropriate experimental investigations of the problem, 
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STABILITY OF SOME CARBONATES AT 25°C 
AND ONE ATMOSPHERE TOTAL PRESSURE* 
RK. M. GARRELS. THOMPSON, and R. SIEVER 


Division of Geological Science, Harvard University, Cambridge, Massachusetts 


ABSTRACT. The stabilities of a number of carbonates, both minerals and synthetic com- 
pounds, were determined by dissolving them in water saturated with one atmosphere car- 
bon dioxide at 25°C and one atmosphere total pressure, The pH of the solution was 
monitored continuously and the equilibrium value determined by extrapolating the pH- 
time data to infinite time. This equilibrium pH value was used to obtain the standard 
free energies of formation of the carbonates, The values obtained are listed below, and are 
compared, where data are available, with values listed in U. S, Bureau of Standards Cir- 
cular 500 or in Latimer’s “Oxidation Potentials.” 


Carbonate This work lie U.S. Bureau 
of Standards 

Calcite, natural and synthetics 269 Be 269.78 
Magnesite 246, 
Dolomite 520.05 
Huntite 1007.7 
Hydromagnesit« 1108.3. 
Rhodochrosite, natural 195. 195.4 
Rhodochrosite, synthetic pptd 194, 194.3 (Latimer) 
Kutnahorite 66.2 
\ragonite 269.6, 269.53 
Witherite 272.0, 272.2 
Alstonite 543.05 
Barvtocalcite 522.9 


Strontianite, natural and synthetic 272.0. 271.9 


The value obtained for dolomite indicates that dolomite is stable with respect to cal- 
cite plus magnesite: conflicting experimental results by previous workers apparently are 
explained by the greatly increased solubility of dolomite that results from prolonged 
grinding 


INTRODUCTION 

One of the oldest methods of determining the solubility of a carbonate is 
to dissolve it in a water solution under a fixed pressure of carbon dioxide, The 
justification for reporting the present work comes from the use of continuous 
recording of pH of the solution. which provides data that can be used to ob- 
tain accurate equilibrium values by appropriate extropolation, Equilibrium 
was approached at 25°C. and one atmosphere pressure only from undersatura- 
tion: but the values we obtained agree well with previous determinations in 
which equilibrium was approached from both sides. If the approach to equi- 
librium from undersaturation alone is valid, the stability of a number of 
species such as dolomite and huntite can be determined, even though these 


compounds have not yet been synthesized at room temperature and pressure. 


Published under the auspices of the Committee on Experimental Geology and Geo- 
physies and the Division of Geological Science at Harvard University. 


Natural or synthetic not specified. 
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The method used has the additional virtue of being rapid and simple; 24 hours 
is usually sufficient to obtain data for determination of an equilibrium constant. 


EXPERIMENTAL 


Vaterials.Water was prepared by distillation in a Barnstead still, and 
yielded a pH of 7.0 + 0.1 after washing with nitrogen, Carbon dioxide was 
obtained from a commercial tank, and was bubbled through a gas-washing 
bottle containing distilled water to presaturate it with water before it entered 
the experimental container. Nitrogen, used for removing other gases from the 
water, was obtained from commercial tanks, and was passed through a carbon 
dioxide absorbent and through distilled water before use. The carbonate com- 
pounds were either mineral specimens of good quality or Reagent Grade com- 
mercial products (table 1). We are indebted to C, Frondel and C, Hurlbut of 
the Department of Mineralogy, Harvard University, for most of the specimens, 
Jun Ito, of the same Department, for analyses of magnesite, huntite, and dolo- 


TasLe | 
Data on carbonate specimens studied 

Specimen Identification * Analysis sample Locality 
Calcite H.M. 87189 Iceland 

Natural 
Calcite** Baker# Lot Same lot 

Synthetic No, 11246 
Dolomite H.M. 105064 Same sample Oberdorf, 
Styria, Austria 
Tea Tree Gully, 
S. Australia 
Crestmore, 
California 
Oberdorf, 
Styria, Austria 
Soda Springs, 
Idaho 


Huntite** H.M. 106304 Same locality 


Huntite** 106589 Same locality 
Magnesite H.M. 105090 Same sample 
Hydromagnesite : None 


Rhodochrosite H.M, 96030 


Natural 
Rhodo« hrosite 
Natural 
Rhodochrosite 
Synthetic 
Rhodochrosite 
Synthetic 
Kutnahorite 
Aragonite 
Alstonite 


Barytoc alcite 
Witherite 


Strontianite 
Natural 
Mrontianite 

Synthetic 


H.M. 89794 
Baker Lot 
No, 90504 


H.M. 85670 
H.M, 84357 
H.M. 76361 
H.M. 80399 
H.M, 88199 
H.M. 93448 


Baker Lot 
No. 8040 


Harvard Museum 


** Less than 2 


microns 


Same locality 
Same specimen 


Same lot 


Made from analytical 


grade reagents 
Same sample 


Same locality 


Same lot 


# Baker Chemical Company, Phillipsburg, New Jersey 


John Reed Mine, 
Alicante, Colo. 
Franklin, N, J. 


Franklin, N. J. 
Bilin, Bohemia 
Alston Moor, 
England 

Alston Moor, 
England 
Northumber- 
land, England 
Ahlen, 
Westphalia 


i 
. 
Fy. 
‘bate 


4A R. M. Garrels, M. E. Thompson, and R. Siever—Stability of 


mite. Donald L. Graf of the Illinois Geological Survey gave us hydromagnesite. 
Analyses of the other carbonates are quoted from the literature; in a few 
instances they represent analyses of the same specimen, but some of the quoted 
analyses are of other specimens from the same locality (table 2). With the 
exceptions of huntite, hydromagnesite, and the synthetic materials, all of the 
samples used were selected from good quality, coarsely crystallized specimens. 
X-ray diffractometer charts were made of all the samples used, and all were 
found to be pure phases, properly labeled, except for the huntite specimen from 
Crestmore, California, which was found to be about half huntite and half 
aragonite. 
TABLE 2 
Analyses of minerals and synthetic materials used 
rhodo hrosite 
dolomite magnesit« Colorado rhodochrosite kutnahorite 


H.M. 105064 H.M. 105090 Wherry & H.M. 89794 H.M. 85670 
anal: J. Ito anal: J. Ito Larsen, 1917 Frondel, 1955 ~—- Frondel, 1955 


CaO 30, 0.10 0.28 2.08 27.44 
McO 21.78 7.2 0.33 0.69 2.21 
MnO aie 05 59.11 58.27 28.31 
FeO 8 1.16 0.27 0.50 
co, 1.31 1.83 37.89 38.71 41.80 


0.82 


100.12 9959 100.02 100.26 


hunt te 

Australia I i huntite alstonite 

Skinne kinner, HM. 106589 Alston, Eng. 
1958 8 anal: J. Ite Kreutz, 1905 


16.0 5 17.60 
33 3.42 4.25 
1.54 
90.4 29.41 


smal] 


Total 99.80 


Synjphietic Synthetic Synthetic 
te Rhodochrosite Strontianite 


Lot Baker Lot Baker Lot 


Bak#: 
14246 #90504 


insoluble in HC] 0.010 
chloride 0.010 0.002 
sulfate 0.003 0.003 
nitrate 0.010 
iron 0.002 0.001 
other heavy metals (as Pb) 0.005 
barium 0.002 
zine 0.05 


Mg and alkali salts (as SO,) 0.05 0.12 0.022 


. 
» 
| 
: insol 0.44 
a} 
CaO 
MeO 
1 
ign. loss 
CO, 
Cl tr N.D. 
Fe, Al 
nsol nil N.D 3.95 adh 
Potal 100.8 97.7 99.58 
: 
a3 
in 
Tad 
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The mineral specimens were ground dry in a motor-driven agate mortar 
hefore they were dissolved; most synthetic carbonates were sufficiently fine as 
obtained to give satisfactorily high solution rates. No attempt was made to 
size the ground carbonates; in general all material used was less than 0.1 mm. 

Apparatus.—Solution of the carbonates was carried out in cylindrical 
pyrex containers of 600 ml capacity. Covers for the containers were made 
from 1/4 inch thick plastic plate, cut to fit, and having holes to fit #3 and 
#1 rubber stoppers. Communication with the solution was established by 
means of a set of pH-measuring electrodes mounted in rubber stoppers and in- 
serted through the holes in the cover. A resistance thermometer for tempera- 
ture compensation, a platinum electrode for solution grounding, and a gas inlet 
tube were also introduced through rubber stoppers firmly seated in the holes 
in the plastic cover. Leads from the electrodes, thermal compensator, and solu- 
tion ground were connected to an amplifier and an automatic recorder, In 
order to use a magnetic stirrer the water bath was set upon the stirring motors 
and the reaction vessels were wrapped in flexible plastic material and sunk to 
the bottom of the bath. Thus a seal became unnecessary. because the plastic 
envelope, which stretched about 12 inches above the top of the container, soon 
became filled with whatever gas was fed in, and provided the desired pressure 
of one atmosphere. Most experiments were run in duplicate; the pH of each 
container was recorded every 18 seconds. The apparatus used is well balanced 
for measurements of pH to within 0.01 pH units; we estimate the accuracy of 
individual measurements as about + 0.02 pH units. The temperature of the 
water bath was controlled at 25 + 0.1°C, 

Procedure._The electrodes were first calibrated in pH 7 and pH 4 stand- 
ard buffer solutions; the check was always within 0.01 pH units. Then the 
pyrex container and electrodes were carefully cleaned and 500 ml of distilled 
water added. A teflon coated magnetic bar was used to stir the solution, Then, 
if the hydrolysis (i.e., carbonate and water with no CO. gas) pH was to be 
measured, the water was washed with N. to a pH of 7.0 + 0.1, Finally several 
grams of fine grained carbonate were added. The container was encased in the 
plastic envelope, weighed down, sunk to the bottom of the water bath, and 
stirred continuously through the glass bottom of the bath. Then CO, was 
bubbled continuously through the solution. The pH, originally high as a result 
of the hydrolysis of the carbonate, dropped rapidly to a low value, then started 
to rise again as the carbonate dissolved. Several runs in the absence of car- 
bonate showed that the water quickly became saturated with CO., and gave 
pH values consistent with equilibrium between gas and solution (pH = 3.91). 


No attempts were made to compensate for variations in atmospheric pres- 
sure, and hence CO, pressure in the solution. Records of barometric pressure 
during the experiments, obtained from the U. S. Weather Bureau, showed that 
the variations from one atmosphere were not sufficient to affect our calculations. 


The pH values obtained during each run were plotted on linear or semi- 
log paper against time”? to obtain an extrapolation to an equilibrium value at 
infinite time. These plots were empirically chosen because they yielded nearly 
linear plots as equilibrium was approached. 
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After a given pair of samples had run long enough to give a satisfactory 
approach to equilibrium, stirring was stopped, and the supernatant solution 
and suspended carbonate removed by decanting. Fresh water was added, and 
the run repeated. This procedure was used to avoid the possibility of a signifi- 
cant contribution to the solubility from material so fine grained as to have a 
significantly higher surface energy than bulk material. 


Most runs required 12 to 18 hours under these experimental conditions. 


When some of the more slowly soluble carbonates required more time, we tried 
to bring them back into the 12 to 18 hour range by using larger samples or 
finer grinding. or both. If the runs lasted 3 or 4 days, contamination of the 
solution and of the calomel electrode resulted, and there was difliculty in ob- 
taining checks on buffer solutions. Also, diffusion of KCI from the calomel 
electrode caused an unknown increase in ionic strength of the solution, 

Analyses of solutions that had approached equilibrium were made from 
time to time to check calculated values, The concentration of bicarbonate ion 
was determined by titrating the filtered solution with standard HCl under one 
atmosphere CO. to an endpoint at pH 3.91, Calcium was determined by precipi- 
tation as the oxalate, or calcium and magnesium were determined by titration 
with Versene, according to the procedure outlined by Shapiro and Brannock 
(Shapiro and Brannock, 1956). Barium was determined gravimetrically as the 
sulfate. and manganese colorimetrically after oxidation to permanganate. The 
results are discussed under the sections devoted to the behavior of the individ- 
ual carbonates. 

RESULTS 

General._-lypical results of the experiments are shown in figure 1. 
Figure la shows curves for huntite and calcite. The similarity in behavior in- 
dicates that a similar solution mechanism is involved, despite the fact that 
huntite contains two different essential cations. Complete equilibrium was 
seldom achieved, but figure 1b, which shows data for synthetic SrCO.,, il- 


Time In Minutes 
2500 400 100 as 25 


T T 


1 L 
05 10 AS 
(Time) 
Log pH versus time for runs of huntite and calcite. 
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Time In Minutes 
1000 400 100 ad 25 16 it 
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Fig. 1. B. pH versus time! for two runs of synthetic strontianite, 


lustrates achievement of equilibrium during the experiment. It would be ex- 
pected, therefore, that our extrapolated equilibrium values from runs of the 
type shown in figure la would give a pH that is too high, and from this a 
lower than actual stability. However, our experimental points are sufliciently 
close to the ordinate to preclude the possibility of a serious error, Also, com- 
parison of our results with those of previous investigators shows that our de- 
viation from them is in the direction of greater, rather than lesser, stability. 

Celcite.-Both natural and synthetic calcites were run. The natural ma- 
terial was large clear rhombic cleavage fragments of Iceland spar, Synthetic 
calcite, obtained from D. Stewart of the U.S. Geological Survey, was a speci- 
men of reagent grade material produced by J. T. Baker Chemical Company. 
Both the natural and synthetic calcite dissolved very rapidly, and a value for 
the equilibrium pH accurate within 0.01 or 0.02 pH units could be taken from 
the highest pH achieved in most of the runs. The filtered solution from some of 
the runs, left open to the room so that the CO. might pass off. precipitated a 
mixture of calcite and aragonite. Equilibrium pH is 6.02. 

Magnesite—The magnesite used was unusually pure material, as shown 
by the analysis in table 2, but a reliable extrapolation to equilibrium pH was 
not obtained because the mineral dissolved so slowly that the curves approached 
the “infinity axis” much too steeply. Numerous runs, using up to 30 grams of 
finely ground magnesite, were made, with inconclusive results. The filtered 
outgassed solutions formed no precipitate for months, until, because of evapora- 
tion, nesquehoneite, MgCO,.5H.O formed. A cleavage fragment of magnesite, 
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suspended in one such solution for 18 days. showed no gain or loss of weight. 
No satisfactory equilibrium pH was obtained. 


Dolomite. The dolomite specimen used was also of very pure material 


(see analysis in table 2). Unlike the magnesite from the same locality, how- 
ever, dolomite behaved very well when subjected to our standard procedure. 
Good duplication of results was obtained, even from the earliest runs, made 
outside the water bath. However. if a run was allowed to approach to within 
about 0.1 pH units of the extrapolated equilibrium value (fig. 2), and then 
subjected to the addition of several grams of dolomite that had been ground 
for 12 hours or more. the pH rose rapidly, and leveled off to give an extra- 
polated value appropriate to aragonite. When this experiment was re-run, after 
Time In Minutes 

178 100 64 44 


10,000 1600 400 


_-—— Addition of 3 grams of dolomite 
a ground for 18 hours. 


025 050 075 100 125 4150 
(Time)"'/2 


Fig. 2. Effect of adding finely ground dolomite to a solution approaching saturation 
with coarser material, 
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removing the supernatant liquid and its accompanying fine particles, an ex- 
trapolation to the usual pH was obtained, Equilibrium pH (coarse material ) 
is 5.08. 

Huntite.Two specimens of huntite were run, but that from Crestmore 
was found to be a mixture of about equal proportions of aragonite and huntite. 
That from Australia was a single phase, however, and good duplication of re- 
sults was obtained from it. Moreover, analyses of some of the solutions showed 
that the mineral dissolved essentially congruently, However, we cannot be sure 
that huntite is as stable as calculated; it is possible that it dissolved congruently 
until calcite started to precipitate, and then ceased reacting. The precipitate 
that formed when some of the filtered solution was left open to the room was 
calcite. 

The runs of the Crestmore aragonite-huntite showed an interesting break 
in slope, from which two extrapolations might be made—the lower one to ap- 
proximately the equilibrium pH of aragonite, and the upper to about the 
equilibrium pH of huntite. Equilibrium pH is 6.24. 

Hydromagnesite-—The hydromagnesite specimen, from near Soda Springs, 
Idaho, was not analyzed chemically, but the x-ray diffractometer chart showed 
it to consist essentially of hydromagnesite, with no other carbonates present, 
but with 10 or 15 percent of silt and clay. Two runs of this material gave 
duplicate results, but the standard HCl titrations of the solutions did not con- 
firm the measured pH’s. It may be that at pH’s of 6 or higher, the Mg* * ion 
is partially complexed, and the measured pH does not give a true measure of 
the cation concentration, Comparison of a series of calculated versus measured 
bicarbonate molalities for magnesium bicarbonate solutions indicate the pres- 
ence of a minor amount of associated Mg(HCO,). or Mg(HCO,)*. We re- 
vard the equilibrium constant value for hydromagnesite as the least reliable 
of the listed values. Equilibrium pH value is 6.71. 

Rhodochrosite. Two mineral specimens of rhodochrosite were run; both 
contained about 5 mole percent of other cations in diadochic substitution. Good 
extrapolations to an equilibrium pH were not easily obtained from either 
rhodochrosite, so the average result of a number of runs is presented. Precipi- 
tated MnCO., prepared from Reagent Grade reagents, was shown by the x-ray 
diffractometer to have the structure of rhodochrosite. though the reflections 
were broad, This material was found to be more soluble than the mineral 
rhodochrosites. As differences in solubility were not found between other min- 
eral and synthetic carbonate pairs, it may be that the difference found for 
MnCO, is due to the presence of about 5 mole percent of other cations in the 
natural rhodochrosites. A small amount of better crystallized synthetic rhodo- 
chrosite was made by using an electric current to cause reaction of Mn* 
and HCO,~ ions in a salt bridge. The rate of precipitation was necessarily 
slow, and only about 100 mg were obtained, (The reaction between Mn* 
and HCO,~ produces CO, as well as MnCOs, but the side effect of oxidation 
of Mn** is less at the lower pH of the HCO,~ solution than at that of the 
higher pH of a CO,~~ solution.) One equilibrium pH run only was made with 
this material because so littke was available, and a projected equilibrium pH 
value of 5.73 was obtained, Equilibrium pH of natural rhodochrosite is 5.09; 
of fine grained precipitated MnCO, is 5.35, 
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AKutnahorite.—Our sample of kutnahorite was part of an analyzed speci- 
men from Franklin, New Jersey. The x-ray diffractometer chart of our sample 
showed only kutnahorite peaks, but the solution behavior suggested that it 
might be contaminated with a little calcite, or perhaps the grinding involved 
in the preparation of the sample tended to disorder the cations, The first few 
runs projected to and reached considerably higher pH’s than later runs, The 
filtered solution from one of the later runs was titrated with standard acid, and 
confirmed the measured pH exactly. The colorimetric analysis showed that Mn 
was 49.5 mole percent of the total cations. Equilibrium pH is 5.58. 

fragonite._The aragonite sample was prepared from beautifully crystal- 
lized material from Bilin. Czechoslovakia. Two runs were made, giving good 
duplication of results, and checking the Bureau of Standards value very close- 
ly. Equilibrium pH is 6.08. 

Witherite—The sample of witherite was prepared from one large crystal. 
Two runs were made, giving good duplication, and titration of the solution for 
HCO and precipitation of Ba as the sulfate both confirmed the measured 
pH. Witherite precipitated again from the filtered solution, Equilibrium pH is 
5.93 

{/stonite. The alstonite sample was prepared by handpicking to separate 
it from the associated calcite and witherite, Most of the fragments of alstonite 
were recognized by their distinctive crystal habit. as well as by distinctive 
luster, The x-ray diffractometer chart showed no extra peaks. Three runs were 
made, giving reasonably good duplication, Titration of one of the solutions for 
HCO confirmed the measured pH. The phase that precipitated out of one 
of the filtered solutions was a barian calcite. Equilibrium pH is 5.93, 

Barytocalcite.-This sample was also prepared by handpicking from the 
associated barite and calcite. The barytocalcite was recognized by a difference 
in luster, The x-ray diffractometer chart showed no extra peaks. Two runs 
were made, giving reasonably good duplication. 

The precipitate that formed in a filtered solution of barytocalcite was a 
barian calcite that was found by analysis to contain 40 mol percent BaCO,. 
Equilibrium pH is 5.95. 

Strontianite.—Both natural and synthetic strontianites were run. and 
good agreement was obtained between the coarsely crystallized mineral and 
the fine grained synthetic precipitate, Titrations for HCO confirmed the 
measured pH’s, and the phase that precipitated from the filtered solution was 
again strontianite. Equilibrium pH is 5.73. 

Summary.—The equilibrium pH values for the carbonates, as well as their 
dissociation constants. are summarized in table 3 and figure 3. In the figure. 
equilibrium constants for all carbonates are shown on the basis of dissociation 


to a single carbonate ion (e.g. k dolomite a Aco). The posi- 


tion of a given carbonate in the table is a criterion of its stability in aqueous 
solution in the absence of « omplexing agents. 


DISCUSSION OF RESULTS 
Standard free energy values.—The equilibrium constants for the ioniza- 
tion of the carbonates, and hence their standard free energies of formation 
from the elements at 25°C and 1 atmosphere total pressure can be obtained 
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Huntite 


Aragonite 


Calcite 


Barytocaicite 
Alistonite Witherite 


Strontianite 
Dolomite 


Kutnahorite 


Rhodochrosite 
(synthetic) 


Rhodochrosite 
(notural) 


| 
8 9 10 
pK 
Fig. 3. Equilibrium pH for solutions saturated in carbonate compounds versus dis- 
sociation constant of solid carbonate, 25°C and 1 atmosphere COs, pressure. 


from the equilibrium pH values. The calculations depend upon AF;° values 


already available in the literature for various ionic and molecular species. so 
that the values obtained for the carbonates are no better than the other data 
used. 


TABLE 3 
Equilibrium pH values for carbonate species 
Equilibrium pH Species 


5.09 rhodochrosite, Colo, 

5.09 rhodochrosite, N. J. 

5.35 svn. MnCOs, 

5.58 kutnahorite 

5.70 dolomite 

5.73 s_rontianite 

5.73 syn, SrCOs 

5.93 witherite 

5.93 alstonite 

5.95 barytocalcite 
calcite 
syn. CaCOs 
aragonite 
huntite, Australia 
hydromagnesite 
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4 sample calculation for dolomite should suffice to show the procedure 
used in obtaining AF,° values. 


For the dissociation of dolomite, CaMgi(CO,). Ca++ + Mgt 
2CO , the equation for the solubility product, K, is: 
The value for aco,-- is obtained from the equilibrium pH through the 


following relationships from the equilibrium for carbonic acid (Harned and 
Owen, 1958). 


Pco an, CO, (2) 
ane Anco an, Co (3) 
Aco moa,” (4) 
Combining equations #2 and #3, and setting Peo, I: 
10-*35 x = (5) 


The values for aca++ and ay,++ are also obtained from the equilibrium pH, 
though less directly. When dolomite dissolves to equilibrium in water saturated 
with CO, at one atmosphere pressure: 

CaMg(CO,;), + 2H,0 + 2CO, = Ca++ + Mg++ + 4HCO,— (6) 

At the pH values (5-7) found experimentally’, no dissolved ionic species 
other than Ca++, Mg++, and HCO,~ need be considered as being of quanti- 
tative importance in the mass balance, hence a close approximation is: 


2Meat+ + Myce, (7) 


Because dolomite apparently dissolves congruently (the ion products of 


Ca++ or Me times CO are less than those for any other known solid 
carbonate spec ies): 
Me» My 4Amuco (8) 


If the ionic strength is known. then ionic activity coefficients (y) may be 
calculated and the value for ayeo-, from equation #5 may be converted to 


Muco, -Muco, , in equation #8, yields values for mea++ and myg+, which 
are in turn converted to acy++ and ayer. The equation for ionic strength is: 
Again, the only quantitatively important species in the solution are the 
dissolved cations and HCO, ~. Therefore, in the case of dolomite: 
1/2 (4mc, + + muco,-) (10) 


In general. for doubly charged cations, the expression for p» is reducible 


1/2 (3muco.” ) (11) 
Using the value for » in the Debye-Hiickel equation, (Klotz, 1950) 
log A z7Vp/(1 + BY (12) 
values for yuco... year+; and are obtained. 


An analytical determination of myeo; may be made by titrating the solu- 
tion with a standard acid, but in general myeo, was calculated from ayco 


? At higher pH values the approximation given by (7) is not so close. 
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by a method of successive approximations, For the first approximation aycos- 
is assumed to be equal to mycos-- From this, ionic strength and yuoos- are 
calculated. Then yucos- is used to correct aycos- and new values for ionic 
strength and yucos- are calculated. In most cases it is necessary to repeat this 
cycling only two or three times. An example of the calculation for dolomite is 
given. To simplify the calculations y’s are read from a previously prepared 
graph of y; versus » taken from a table in Klotz (p. 332). 
Sample calculation.—E-quilibrium pH for dolomite = 5.68 (table 3). 

ancos- = = 0.00724 
Muco aHcos 0.00724 (First approximation ) 
= 3,2 (0.00724) = 0.01084 

yuoos 0.902 
Myce ancos-/ 902 0.00724, .902 0.00802 (Second approximation ) 

= 3/2 (0.00802) = 0.01202 

yucos- = 0.898 
0.00724/.898 0.00805 (Third approximation) 

pw = 3/2 (0.00805) 0.01208 

yucos = 0.898 

1/4 mycos- 0.00201 

Yeas = 0.654 

YMg** 0.665 
Acar (0.00201) .654. 0.00131 =-10—?-§§ 
ayg++ (0.00201) .665 = 0.00135 = 


10 5,68 
K (10-2-8) (10—2-87) (10—®79) 2 
From the relation AF,” —RT In K = -1.364 log k at 25°C and 1 at- 
mosphere pressure, the standard free energy of reaction 9 is: 
AF,’ —1.364 X —19.33 26.37 kilocalories (13 
Then we can write: 
Substituting the value for AF,” obtained in equation 13, and values for 
the ions from U.S. Bureau of Standards Circular 500: 
+ 26.37 = —132.18 —108.99 —252.44 —AF 
AF = —519.9, kilocalories 
A check on the consistency of the free energy values from the literature 
can be obtained by calculating AF;° dolomite from the initial reaction with 
CO.: 
CaMg(CO,). + 2H.O0 + 2CO, = Ca++ + Mg++ + 4HCO,— 
kK 10-14-31 
AFR® = +19.52 kcal. 
+19.52 - AF T AF 4AF ’ucos- A camg«cos)2 
—2AF — 2AF 
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19.52 132.18 108.99 — 561.24 —AF;’camgicosy2 + 113.4 + 188.52 
AF cos): 520.0, kilocalories 

Our estimate of the accuracy of the free energy values, based on the 
checks against known values where available, the consistency of values obtained 
by various calculation routes, the generally good agreement between analyzed 
solutions and values calculated for them and the probable error in determining 
the extrapolated equilibrium pH value, is that they are generally accurate to 
about + 0.2 kilocalorie. This assumes no error in the Bureau of Standards 
values for various species used in calculating AF,” of a given carbonate from 
the standard free energy of the dissociation reaction. 

Relations among the Ca-Mg carbonates.—Relations among the Ca-Mg 
carbonates, although not entirely definitive because of our inability to get a 
satisfactory free energy value for magnesite, indicate that calcite and magne- 
site are unstable with respect to dolomite at 1 atmosphere total pressure and 
25°C. Using the U.S. Bureau of Standards value for magnesite: 

CaCO. + MeCO CaMeg(CO,).; AFR” = kilocalories 

Huntite apparently is not stable with respect to dolomite and magnesite: 
CaMe(CO,). + 2MeCO CaMeg,(CO,),; AFR’ + 5.1 kilocalories 

Nesquehonite is markedly unstable (Latimer, 1952). Even in the presence 
of pure water: 

MegCO,.3 H.O MeCO, + 3 AF,’ —4.5 kilocalories 

Thus in the composition triangle CaCO,-MgCO,—H,O, relations among 
the species studied are shown in figure 4. 

The solubility of the carbonates in water is so small that the composition 
of the solutions in equilibrium with the solids cannot be shown well on figure 
t. A more graphic illustration of the relations in this system is obtained by 
plotting stability of the solids as a function of the activities of the Ca** and 


Mg++ ions in solution, and as a function of Peos. This permits inclusion of 
hydromagnesite, which falls outside the triangle CaCO,-MgCO,—H,0. 
All the reactions of interest can be written in terms of the three variables 


Acar, Auge, and Peo, always assuming water in excess, For example, for 
equilibrium between calcite and dolomite: 
2 CaCO, + Mgt + CaMgiCO,). + Cat+ (14) 


For equilibrium between calcite and hydromagnesite: 
CaCO, + 4Me++ + = + 
K 
Me** 

Consequently a logarithmic plot of the variables provides linear bound- 
aries among species. 

The summary diagram is shown in figure 5, Note that hydromagnesite is 
not stable under earth surface conditions (Peo. 10~-*-°), Moreover, if the 
ratios of the molalities of Ca++ and Mg** in sea water are used as criteria 
of their activities, the stable carbonate in sea water should be dolomite. In 
other words, the implication is that the oceans are supersaturated with respect 
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Nesquehonite 


Coco, Comg(CO5), [camg,(CO,), 
Calcite Dolomite Huntite 

Aragonite 


Fig. 4. Some stability relations in the system CaCO;-MgCO,;-H.0. 


Magnesite 


to dolomite, and calcite or aragonite on the floor of the sea should react with 


the sea water to form dolomite until the ratio of the activity of Ca++ to Mg++ 
is equal to 10**-*, 


There appear to be two alternatives to this conclusion, First, our value for 
the free energy of formation of dolomite may be in error, Second, the concen- 
tration of magnesium ion in sea water may be vastly different from its activity. 
The approximate activity of calcium ion in sea water is not in serious ques- 
tion,* therefore any large anomalies between concentrations and activities must 
be attributed to the magnesium ion, The small amount of indirect evidence 
available now suggests that the value we have obtained for the free energy of 
formation of dolomite is approximately correct, and that the near-absence of 
dolomite as a primary marine precipitate results from supersaturation. 

Robie (personal communication) has summarized all the data available 
for calcite, magnesite, and dolomite, and has determined values for their free 
energies of formation that are in essential agreement with ours. 

The solubility data of Yanat’eva (1954) on calcite, magnesite, and dolo- 
mite, as determined at 25°C and 1 atmosphere Peo., show dolomite as a stable 
phase. On the other hand her data at 25°C, 1 atmosphere total pressure but 
with Peo: = 107% (that of the present atmosphere), indicate that dolomite 
dissolves incongruently to yield calcite plus solution. As shown in figure 5 or 
in equation number 14, the stability of dolomite relative to calcite is not a 
function of Peas. so the inference is that her experiments at lower Peo: and 
hence higher pH, somehow caused dolomite to increase its relative solubility 
drastically. 

* The activity coefficient for Ca** in sea water is of the order of 0.24, as shown by various 
studies of the solubility of calcite. (c.f. Garrels and Dreyer, 1952, p. 334.) 
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Fig. 5. Stability of some Ca-Mg carbonates in aqueous solution at 25°C, and one 
atmosphere total pressure, as a function of calcium-magnesium ion ratio and partial pres- 
sure of carbon dioxide 


One possible explanation of this peculiar phenomenon is that the behavior 
of dolomite is analogous to that of quartz. When quartz solubility is deter- 
mined by experiments during which the quartz is tumbled or otherwise ground, 
the solubility increases by 10 to 15 times and approaches that of amorphous 
silica (Van Lier, 1959). This increased solubility is considered to result from 
the disordering of the quartz surface during grinding, so that a layer of 
amorphous silica is produced, Because dolomite has two structurally non- 
equivalent cations it is susceptible to disordering effects from grinding. This 
disordering has been reported by Bradley, Burst, and Graf (1953), These ef- 
fects would be especially noticeable in experiments in which the dolomite 
would be continuously stirred or tumbled for long periods of time, The effects 
of a very slight disordering produced during the preparation of the sample 
would be noticed especially in experiments where the solubility is small, as 
under low CO, partial pressures. 

To test this hypothesis we permitted dolomite to react with water saturated 
with CO, under one atmosphere pressure, until the pH of the solution was with- 
in 0.1 pH unit of the extrapolated equilibrium value. (Note that under these 
conditions dolomite is soluble enough so that if any small portion of it has 
been disordered by the original preparation of the sample it will all be dis- 
solved before the solution is saturated.) Then three grams of the same dolo- 
mite, that had been ground continuously for 18 hours, were added to the 
system. The effect was immediate and striking. Within minutes the pH had 
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risen nearly 0.5 pH units, and after a few hundred minutes a new extrapolated 
equilibrium pH was obtained (fig. 2). Calculations based on the higher 
value (pH 6.2) show that the product aca++ * acos-- is 10~***, exactly the 
value we had previously determined as the equilibrium constant for aragonite. 
Apparently the prolonged grinding had so increased dolomite solubility that it 
dissolved incongruently to yield CaCO, plus magnesium-rich solution. 

It was also found that when the first solution was decanted and replaced 
vith more distilled water and the solution experiment repeated that the dolo- 
mite had reverted to its original solubility. In other words, the material re- 
sponsible for the increased solubility had all been dissolved during the first run 
oc was carried off by the decanting, and by the time the solution approached 
saturation the second time, we were again measuring the solubility of relative- 
iy coarse-grained, ordered dolomite. 

The anomalously high solubility was measured again by addition of more 
of the well-ground dolomite to the new solution, with the same results. 

It is concluded that our usual procedure measures the stability of dolomite 
accurately, and that the higher solubilities reported by other investigators re- 
silt from grinding, either during sample preparation or from the stirring or 
tumbling during the measurement itself. 

Our experiments are too preliminary in nature to explain this phenomenon 
completely. One possibility is that the surfaces of the grains are disordered 
during prolonged grinding; another possibility is that a small amount of the 
sample is ground so fine that its solubility is enhanced, Ordinarily, when a 
supersaturated solution is in contact with the appropriate solid, the supersa- 
turation is relieved by precipitation on the less soluble larger grains, but be- 
cause of the extremely slow rate of precipitation of dolomite, this effect is 
probably not seen in the usual experiments, X-ray diffractometer charts of the 
dolomite before and after grinding showed that the prolonged grinding tended 
to produce line broadening due to decreased particle size, There was no ap- 
parent relative intensity decrease of the ordering reflections. 

Similarly, well-ground calcite was added to a calcite run near the equilib- 
rium pH. After the first addition the solution attained the equilibrium pH and 
successive additions of finely ground material raised the pH by no more than 


0.01 pH units. 


Within 5 minutes of each addition, the pH was again at the equilibrium 
value. The effect in this case is so slight that we cannot definitely say whether 
the small increase in pH was due to excess solubility of the finely ground ma- 
terial, or whether it was due to a slight loss of CO, atmosphere as a result of 
the necessary opening of the container. 


Furthermore, our whole experience with magnesium ion in solution has 
impressed us with its reluctance to form solids irom saturated solution, When 
solutions partly saturated with magnesite at Peo: = 1 atmosphere were allowed 
to stand open, so that the Peo: approached that of the laboratory, no precipi- 
tate formed over a period of several months, These solutions were highly super- 
saturated with respect to magnesite. If precipitation was finally forced by 
permitting the solutions to evaporate, the product was nesquehonite, which is 
many times as soluble as magnesite. 
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Perhaps the clue to the reluctance of magnesian carbonates to precipitate 


lies in some peculiar characteristic of the magnesium ion; it may be abnormal- 
ly firmly hydrated, or it may form a stable ion pair (MgCO,)°. 
Calcium-barium carbonates.—-The data on the calcium-barium carbonate 
system indicates that the ordered intermediate phases barytocalcite and alston- 
ite have approximately the same solubility, and that they are slightly more 
stable than an ideal solid solution of the same composition. Our data are not 


good enough to permit us to determine which of the dimorphs is more stable. 
but the free energy of the reaction 


BaCa(CO,). (Barytocalcite ) BaCaiCO,). (Alstonite) 


is probably not more than + 100 calories. In harmony with our other attempts 
to form ordered intermediate phases, the precipitate resulting from loss of CO. 
from a nearly saturated CO,-rich solution yielded a high (40 mol %) barian 
calcite, rather than the initial barytocalcite or alstonite. 


Ca-Mn_ carbonates.Results of experiments involving kutnahorite indi- 
cate that this ordered intermediate phase is stable at the conditions of the ex- 
periment. Also. we observed a marked difference in the stability of precipitated 
and natural MnCO,. The difference we found is about the same as that re- 
ported by Latimer (Latimer, 1952, p, 235) 


Sr.-carbonate.—Our value for strontianite agrees closely with that of the 
Bureau of Standards, No anomalies in behavior were observed. 
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CLASSIFICATION OF ARENITES 
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University of New England, Armidale, N.S.W., Australia 


ABSTRACT. A critical examination of the classificatory schemes of Pettijohn (1954), 
Gilbert (1955), Packham (1954), and Folk (1954) leads to the conclusion that only 
Packham’s, which uses sedimentary structures as a parameter, is validly based genetically. 
Packham’s scheme is expanded by the introduction of a generalized, non-genetic, ‘arenite’ 
diagram which is congruent with the genetic graywacke and arkose-quartzose sandstone 
(here termed simply ‘sandstone’) diagrams and which enables it to deal with arenites 
when data on sedimentary structures are not available. 

A new QFR diagram is introduced having quartz, feldspar, and rock fragments plus 
other labiles as parameters, The non-genetic classes of labile arenites recognized are 
feldspathic, lithofeldspathic, feldspatholithic, and lithic. Congruent graywacke and sand- 
stone QFR diagrams may be employed where appropriate. Two classes of sublabile arenite 
are recognized—feldspathie sublabile arenite and lithic sublabile arenite, The term 
‘arenite’ may again be replaced by ‘sandstone’ or ‘graywacke’ where appropriate, The use 
of ‘feldspathic sandstone/arenite/wacke’ as descriptive of sublabile arenites should be 
discontinued, 


INTRODUCTION 
Since the publication of Krynine’s classic paper on sedimentary rock 
classification (Krynine, 1948), in which ternary diagrams were first utilized 
for classifying arenites, a number of classificatory schemes for detrital arenites 
have appeared. Those of Pettijohn (1949), replaced by his scheme of 1954, 
Gilbert (1955). Packham (1954). and Folk (1954) have found the widest ac- 
ceptance, and have largely superseded Krynine’s scheme, Certain shortcomings 
exist in each of the newer classifications, but, as will be shown, only those in 
Packham’s can be remedied without destroying the classification. 
Concerning the aim of classifications, Pettijohn (1948) writes: 
No classification independent of rock genesis can be worthy of consideration. 
Genesis is the ultimate aim of any study of rocks, and no descriptive classifica- 
tion—so-called—can be worth much unless the characters used for the classi- 


fication are meaningful and significant. The only test for significance is whether 
they are or are not basic to the understanding of origin. 


It would seem proper, therefore, to judge the current classifications on the ex- 
tent to which they are truly genetic. 

A useful way of assessing the genetic basis of a given classification is to 
classify a group of arenites deposited under similar conditions, and compare 
the names given with concepts of genesis derived from other sources, For ex- 
ample, graywacke and orthoquartzite are commonly held to be characteristic of 
different tectonic and environmental settings. The occurrence of these two rock- 
types in one formation would suggest most unusual conditions or a certain 


arbitrariness in the classificatory scheme used. If, after several analyses, a 
continuous gradation were found between graywacke and orthoquartzite, the 
latter alternative could be asserted with confidence. 


TERMINOLOGY 
Arenite.-Arenite is used in the literature in two senses. The most com- 
monly accepted is that of Pettijohn (1957, p. 16) ; a descriptive term for clastic 
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material of sand size without genetic or mineralogical connotations. All clastic 
rocks composed of sand-sized fragments are arenites in this definition. 
“arenite’” to a group 
of sand-sized sediments which contain less than 10 percent detrital argillaceous 
material, Gilbert considers “‘arenites” to be genetically distinct from ‘wackes’ 
which contain more than 10 percent detrital argillaceous material. 

In the proposed expansion of Packham’s classification which follows, the 
term ‘arenite’ is used in the broader, non-genetic, sense of Pettijohn, This 
usage is followed throughout the paper, except when discussing Gilbert's 
“arenites”. In such cases the change in sense is apparent from the context. 

VLQ and QFR diagrams.—Two types of triangular diagrams have been 
widely employed in the classification of arenites, Pettijohn (1949) and 
Packham (1954) have used MLQ diagrams which have as parameters per- 


Gilbert (1955, p. 289) proposed limiting the term 


centages (on a cement-free basis) of matrix (M), labiles (feldspar, rock frag- 
ments, ferromagnesians, chlorite) (L), and quartz plus chert (Q). 

Gilbert (1955) on the other hand, utilized QFR diagrams, whose para- 
meters are percentages (on a matrix- and cement-free basis) of quartz plus 
chert (Q), feldspar (F), and rock fragments plus other labiles (R). 

Detrital carbonate is excluded from both types of diagram. Arenites con- 
taining abundant detrital carbonate require other treatment. 


DISCUSSION OF CLASSIFICATIONS 


Pettijohn’s 1954 classification —Pettijohn (1954) recognized three factors 
to be of greatest genetic significance: a maturity factor (ratio of quartz plus 
chert to labiles on a matrix-free basis) ; a provenance factor (ratio of feldspar 
to rock fragments), which he considers to be a measure of the plutonic and 
supracrustal rocks exposed in the source area; and a fluidity factor (ratio of 
sand detritus to martix), which gives an indication of the viscosity of the de- 
positing medium. His classification, which is shown in the form of an MLQ 
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Fig. 1. Comparison between Pettijohn’s 1949 and 1954 classifications. 
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diagram in figure 1B, is based on these factors. A noteworthy feature is the 
location of the graywacke field boundary at a 15 percent level of matrix. 

Criticisms.—First, as originally defined, Pettijohn’s provenance factor 
measures the nature of the source area only in the special case where all felds- 
par is plutonic and all rock fragments are supracrustal, If this factor is applied 
to the Miocene graywackes of the Aure Trough, Papua, described by Edwards 
(1950). for example, it leads to completely erroneous conclusions, for these 
rocks are dominantly of andesitic derivation. This difficulty seems to have 
been realized subsequently by Pettijohn (1957, p. 290), as his phraseology is 
somewhat modified. 

Second, the choice of the boundary of the graywacke field, while not with- 
out reason, is arbitrary. The arenites of the Burralow Formation, Triassic, of 
New South Wales (Crook, 1957, p. 69), for example, belong to a formation 
that was deposited under broadly similar conditions, but they lie over the 
graywacke, protoquartzite and orthoquartzite fields of Pettijohn’s diagram. 

Packham (1954, fig. 1) has plotted the frequency of matrix percentages 
for arenites deposited from traction and turbidity currents, An overlap is 
shown between 8 and 32 percent of matrix, so that Pettijohn’s decision to 
characterize graywackes by a minimum level of 15 percent matrix must be 
recognized as a good first approximation, but nonetheless arbitrary. Its use 
destroys the truly genetic nature of the classification, since the occurrence of 
cases such as the Triassic example cited is inevitable. This point has been well 
made by Folk (1954, p. 353), in discussing the confusion between texture and 
composition in certain classifications: “ . . , one’s terminology would fluctuate 
back and forth . . . merely at the whim of the depositional agent.” 

A final criticism lies in the changed definition of “subgraywacke’ from 
the 1949 to the 1954 classifications (see figs. LA and 1B). This change, while 
bringing the term more into line with Continental usage, has added to the con- 
fusion in nomenclature. 

Gilbert's classification.—Gilbert’s classification is basically similar to 
Pettijohn’s 1954 classification; it employs the same three factors, by implica- 
tion, but uses two diagrams of the QFR type. The first (fig. 2A) is the ‘wacke 
diagram’ for arenites with more than 10 percent detrital matrix, Thus it is the 
equivalent of those parts of Pettijohn’s 1954 MLQ diagram (see fig. 1B) which 
lie above the 10 percent level of matrix. The second QFR diagram (fig, 2B) is 
the ‘arenite diagram’, for arenites with less than 10 percent detrital matrix, It 
is the equivalent of the remainder of figure 1B. 

Criticisms.—Gilbert’s classification can be criticized in one aspect on the 
same grounds as Pettijohn’s 1954 classification. The choice between the two 
diagrams depends on the arbitrary level of 10 percent matrix, Like the 15 per- 
cent level of Pettijohn it fails, and for the same reasons. 

A further criticism lies in Gilbert’s usage of the term ‘arenite’, It is most 
unfortunate that this term, previously without mineralogical or genetic con- 
notations, should have been so restricted, as it destroys the only existing satis- 
factory term for the general description of sand-sized sediments, Pettijohn 
(1957, p. 17) has already commented on the inadequacy of the term ‘sand- 
stone’ for such general descriptions. 
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Fig. 2. QFR diagrams of Gilbert's classification (1954). 


Packham’ s classification. Packham (1954) utilizes, by implication, two 
of Pettijohn’s factors: the fluidity factor and maturity factor, A new parameter 
is introduced, however, sedimentary structures (Packham, 1954, table 1). On 
the basis of sedimentary structures arenites are assigned to one of two groups, 
the ‘graywacke suite’ and the ‘arkose-quartzose sandstone suite’, Each group 


has its own MLQ diagram on which arenites may be plotted in the normal 
manner. The name so derived will be at once both descriptive and genetic. 
Whether Packham’s classification is truly genetic depends on the value of 
sedimentary structures as indicators of the type of depositing current, Data 
obtained by numerous workers during the last ten years strongly suggest that 
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Fig. 3. MLQ diagrams for arenites on Packham’s classification (1954). 
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sedimentary structures can be validly used in the manner Packham has sug- 
gested, 

By using two MLQ diagrams, Packham has avoided the arbitrariness 
found in Pettijohn’s and Gilbert's classifications, and it is impossible for speci- 
mens from a single genetically homogeneous formation to have names implying 
different geneses because of slight variations in the detrital matrix content. 

Criticisms.—Packham’s classification, in its present form, is usable only 
when data on sedimentary structures are available. No provision has been 
made for field and laboratory nomenclature which does not imply genesis. 

The term ‘sandstone’ is used by Packham in two senses, broadly, as a 
synonym of ‘arenite’ as used herein, and in a restricted sense, always with a 
qualifying adjective, to refer to arenites deposited by traction currents, for 
example, ‘labile sandstone’, “quartzose sandstone’. This homonymy causes dif- 
ficulties if one is dealing with arenites without knowledge of their associated 
sedimentary structures. They cannot be descriptively named on Packham’s 
scheme without prejudging the issue, An arenite rich in quartz, but of un- 
known associations, cannot be called a ‘quartzose sandstone’ without implying 
the association of certain sedimentary structures, which are however unknown. 

The classification is inadequate when attempting to name arenites which 
contain large percentages of labile material. For rocks derived from non- 
granitic sources only the terms ‘labile graywacke’ and ‘labile sandstone’ are 
available. 

Folk’s Classification. olk’s classification does not utilize either QFR or 


MLQ diagrams, but uses a ternary diagram having quartz plus chert: feldspar 

plus all igneous rock fragments; and mica plus metamorphic rock fragments 
plus metaquartzite as parameters. The diagram is subdivided in the manner 
shown in figure 4. Unlike the classifications already discussed, Folk’s scheme 
does not utilize matrix-content as a parameter, and the names given to arenites 

in his scheme are thus uninfluenced by the matrix-content of the rock, Al- 
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Fig. 4. Triangular diagram for Folk’s classification (1954). 
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though matrix-content effects the plotting of arenites on MLQ diagrams in 
Packham’s scheme (see fig. 3), the end result is the same, for variation in 
matrix-content does not result in changes in arenite nomenclature. 

Folk’s scheme does not attempt to convey information concerning mode of 
sedimentation, the parameters being chosen with the intent to classify arenites 
according to provenance. The three parameters given above are held to meas- 
ure (a) detritus of sedimentary origin, or detritus which has undergone a 
prolonged period of chemical and or physical modification; (b) detritus of 
igneous origin; and (c) detritus of metamorphic origin. 

Criticisms.—-The most important criticisms of Folk’s classification find 
their source in his choice of parameters. The quartz plus chert parameter is 
not a completely satisfactory measure of detritus of sedimentary origin, because 
it may incorporate an unknown quantity of first cycle plutonic or volcanic 
quartz, and volcanic “chert”, which may not have undergone prolonged 
chemical or physical weathering. Furthermore sedimentary rock fragments 
find no place in his parameters, and their insertion poses great difficulties, 

The grouping of all igneous detritus into one parameter is unfortunate, 
for the mineralogical and tectonic significance of plutonic as against volcanic 
components in a sediment is considerable, and they should be separated in ar- 
riving at a rock name. 

The lack of emphasis on mode of sedimentation is, in the author’s opinion, 
a fundamental weakness of the scheme. It leaves no possibility of distinguish- 
ing an arenite deposited by a turbidity current in deep water from one de- 
posited by fluvial agencies on land in those cases, which are not uncommon, 
where both have a large igneous component. 

The extent of the difficulties raised by the choice of parameters and lack 
of emphasis on mode of sedimentation may be seen by comparing arenites 
derived from (a) an andesitic island arc; (b) a granite mass of low relief be- 
ing covered by a transgressive sea; and (c) a complex orogenic belt with 
prominent exposures of volcanics, It is not unlikely that the arenites derived 
from each source will fall in the “arkose” field on Folk’s diagram (fig. 4). 
Although this contingency can be met by the use of suitable qualifiers as Folk 
has suggested (1954, p. 357), it is rather unwieldy and confusing. 


EXPANSION OF PACKHAM S GLASSIFICATION 


From the foregoing Packham’s classification seems to be the only existing 


one soundly based, genetically speaking. It is proposed to use this classification 
as a basis for an expanded classification which will not be subject to the 
criticisms mentioned above. 

To meet the difficulties encountered in dealing with arenites when data 
on sedimentary structures are not available, the author proposes a non-genetic 
‘arenite MLQ diagram’. This diagram (fig. 5) employs the same parameters 
as the graywacke and arkose-quartzose sandstone suite diagrams of Packham, 
and is subdivided in the same manner on the basis of the maturity factor, into 
‘labile arenite’, ‘sublabile arenite’, and ‘quartzose arenite’. The relationship be- 
tween this new diagram and those of Packham is shown in figure 6. Arenites 
of unknown associations, provided they are not rich in detrital carbonate, can 
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Fig. 5. Generalized MLQ diagram for arenites. 


he named by reference to this diagram, and the names derived will be descrip- 
tive and non-genetic. 

‘Sandstone’, will be restricted to sediments of inferred traction-current 
origin, and ‘arenite’ used as the generalized non-genetic term for sand-sized 
clastic material. This avoids the homonymy in Packham’s usage and his 
‘arkose-quartzose sandstone suite’ can then be termed simply the ‘sandstone 
suite’. 

Amplification of the terms descriptive of labile rock-types is desirable, 
particularly to enable distinction between labile arenites rich in feldspar and 
those rich in rock fragments. This can be effected by using a QFR diagram 
with suitable subdivisions. For simplicity only one, generalized, QFR diagram 
is proposed, which gives a series of names for arenites of various compositions: 
congruent ‘graywacke’ and ‘sandstone’ QFR diagrams may be used if data on 
sedimentary structures are available. 

The QFR diagram proposed (fig. 7) is subdivided congruently with the 
arenite-MLQ diagram (fig. 6) into quartzose arenite and sublabile arenite at 
levels of 90 and 75 percent quartz, respectively. The sublabile arenite field is 
subdivided as discussed below. 

The field containing rocks with less than 75 percent quartz is subdivided 
on the basis of ratio of feldspar to rock fragments + other labiles (the FR 
ratio), this ratio being calculated on a quartz-free basis. The subdivided fields 
are designated as follows, and have the following range of values for the FR 
ratio: 


feldspathic arenite FR—2 to 3 7 feldspathic 


lithofeldspathic arenite FR—3 to 1 f labile arenites 
feldspatholithic arenite FR—1 to 1/3 lithic 


lithic arenite FR—1/3 to 0 f labile arenites 
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The author has followed Pettijohn (1954) in using the term ‘lithic’ to desig- 
nate arenites with abundant rock fragments. 

Labile arenites rich in feldspar should logically be termed ‘feldspathic 
arenites’, Unfortunately the adjective ‘feldspathic’ has been commonly applied 
to certain kinds of sublabile arenites. Packham (1954, p. 475) describes the 
sublabile equivalent of an arkose as ‘feldspathic sandstone’, Gilbert (1955) 
uses the terms ‘feldspathic arenite’ and ‘feldspathic wacke’ for rocks which 
would mostly plot in the ‘sublabile arenite’ field of figure 6. 

It is here proposed that the term ‘feldspathic arenite’ be reserved for labile 
arenites with less than 75 percent quartz and an FR ratio exceeding 3 (see 
fig. 7). The genetically flavored congruent terms, ‘feldspathic sandstone’ and 
‘feldspathic graywacke’, should be reserved for arenites of this composition 
which are members of the sandstone and graywacke suites, respectively. Each 
of these three feldspathic rock-types includes, as derivatives of granitic ter- 
rains, the varieties ‘arkosic arenite’, ‘arkose’, and ‘arkosic graywacke’ respec- 
tively. 

The use of the terms ‘feldspathic arenite, sandstone or wacke’ as descrip- 
tive of any variety of sublabile arenite should be discontinued, Pettijohn, 
realizing the desirability of this, has coined the term ‘subarkose’ which is 
roughly synonymous with his 1949 ‘feldspathic quartzite’ (see figs. 1A and 
1B). The genetic connotation of ‘arkose’ limits the value of this term, and it 
should be restricted to sublabile derivatives of granitic terrains, deposited by 
traction currents. 

To overcome this difficulty it is proposed that two broad classes of sub- 
labile arenite be recognized, ‘feldspathic sublabile arenite’ and ‘lithic sublabile 
arenite’, having FR ratios exceeding and less than 1, respectively (see fig. 7). 
As before, the term ‘arenite’ can be replaced by ‘sandstone’ or “graywacke’ 
when data on depositional environment are available. Subarkose then becomes 
a variety of feldspathic sublabile sandstone, and is, in this usage, synonymous 
with Packham’s ‘feldspathic sandstone’. 


CONCLUSION 


So many classificatory schemes for arenites have been proposed in recent 
years that it seems necessary to justify a further increase in their number. The 
following arguments are cited in defense of the scheme presented here: 

a. It is not basically a new classification, but an expansion of an existing 
classification which removes certain shortcomings present in that classification. 
It should therefore be used in conjunction with this existing classification of 
Packham (1954). 

b. Taken in conjunction with Packham’s classification the present scheme 
enables the naming of an arenite in the field, in the laboratory, and, after 
examination of evidence from sedimentary structures, the assigning of a final 
name with both descriptive and genetic connotations. 


It should be noted that the scheme proposed herein does not attempt to 
assign names to all kinds of arenites. This must await further work. Cases 
which are not covered are, for example, pyroclastic arenites, glacial arenites, 
and arenites with abundant detrital carbonate. 
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EVIDENCE FOR AN ABRUPT CHANGE IN 
CLIMATE CLOSE TO 11,000 YEARS AGO* 


WALLACE S. BROECKER, MAURICE EWING and BRUCE C. HEEZEN 


Lamont Geological Observatory, Columbia University, Palisades, New York 


ABSTRACT. Evidence from a number of geographically isolated systems suggests that 
the warming of world-wide climate which occurred at the close of Wisconsin glacial times 
was extremely abrupt. Surface ocean temperatures in the Atlantic Ocean increased by 
several degrees centigrade and deep-sea sedimentation rates sharply decreased in the equa- 
torial Atlantic. The bottom waters of the Cariaco Trench, a restricted basin in the Carib- 
bean, abruptly stagnated, The pluvial lakes in the Great Basin rapidly shrank from their 
maximum volume to nearly their present size, Silty clay which was flushed into the Gulf 
of Mexico during the glacial period was suddenly retained in the alluvial valley and delta 
of the Mississippi River and the glacially-dammed Great Lakes drainage was returned to 
the northern outlet, As indicated by pollen profiles a pronounced warm period in north- 
western Europe also occurred at this time. In each case the radiocarbon age determinations 
suggest that the changes occurred in less than 1000 years close to 11,000 years ago. 


INTRODUCTION 

The problem of reconstructing the variations in climate which characterize 
the Pleistocene has long challenged the geologist. Although the collection, cor- 
relation and interpretation of data has progressed steadily during the past 100 
years, the recent development of geochemical methods for both time and tem- 
perature measurement have greatly accelerated the progress in this field, The 
C** method of age determination in particular has provided not only absolute 
estimates of time over the past 40,000 years but has also allowed precise cor- 
relations to be made between events in widely separated geographic areas and 
in vastly different depositional environments. 

This paper summarizes the evidence for an abrupt world-wide change in 
climate close to 11,000 years ago marking the end of the Wisconsin glacial 
period, That this climatic change occurred suddenly was first recognized by 
Ewing as a result of the study of deep-sea cores (Ewing and Donn, 1956). 
Subsequently radiocarbon measurements on these cores were presented by 
Ericson, et al. (1956). The study has since been extended to other climate- 
controlled systems. 

The evidence to be summarized is largely correlated on the basis of radio- 
carbon measurements. These include not only those reported in “Lamont 
Natural Radiocarbon Measurements” (Broecker, et al., 1956, and Broecker and 
Kulp, 1957) but also the results published by other laboratories. Since the 
period of interest spans less than 2000 years, careful attention must be given 
not only to the laboratory errors, which range between 100 and 500 years, but 
also to the possibility of systematic errors for samples of different material 
types and growth environments. An attempt will be made to evaluate these 
factors as the discussion proceeds. 

The systems studied include deep-sea sediments from the Atlantic Ocean 
and adjacent seas, deposits from the pluvial lake area of the western U, S, A., 
sediments from the Great Lakes and their associated drainage .etworks and 


* Contribution no. 402. Lamont Geological Observatory, Columbia University, Palisades, 
New York. 
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the pollen sequences of northwestern Europe. Each of these systems shows a 
strong dependence on climate. 


EVIDENCE FROM DEEP-SEA SEDIMENTS 

Surjace temperatures. The first and most important of these systems is 
the ocean, Its vast extent and fairly rapid internal mixing rate enable it to 
record world-wide climate changes, being less influenced by local events than 
any other system, This broad picture is recorded in a number of different ways. 
Three of these have been investigated in the course of this study: 1) variations 
in surface ocean temperature, 2) changes in the rate of deep-sea sedimentation, 
}) stagnation of restricted basins. 

A record of surface ocean temperature is provided by the tests of plank- 
tonic foraminifera which make up the majority of the coarse fraction (>74y) 
of the carbonate material in globigerina ooze. This record can be read by either 
the micropaleontological method used by Ericson and Wollin (1956) or the 
oxygen-isotope method used by Emiliani (1955). The first method is based on 
the temperature dependence of the relative abundance of certain species of 
planktonic foraminifera (Schott, 1935). By observation of the relative abun- 
dance of these species an estimate of the relative surface water temperature for 
the time corresponding to any depth in the core can be made. The second 
method relies on the temperature dependence of the oxy gen-isotope fractiona- 
tion that occurs during the formation of carbonate, The isotopic composition 
of the oxygen from carbonate samples can be measured with a precision 
equivalent to 1°C (Urey, 1951). The uncertainty in the absolute temperature 
computed from these ratios is somewhat larger than 1°C because the oxygen- 
isotope ratio for the sea water from which the carbonate was precipitated must 
be estimated from other considerations. For example, Emiliani (1955) has 
pointed out that the preferential storage of O'° in glacial ice could have caused 


a significant change in the ratio for seawater during Wisconsin time (equiva- 


lent to at least 2°C change in temperature}. 

As shown in figure 1 the two methods have shown excellent agreement 
where they have both been applied to the same core. Both methods show a 
rather sharp increase from the relatively cool temperatures of Wisconsin time 
to the warm temperatures of the present. The magnitude of this change as 
given by the paleotemperature method is 6° to 10°C (Emiliani, 1955), de- 
pending on the assumptions made about the glacial storage effect. 

Since the span of time required for the temperature change is critical to 
the arguments which follow, it is important to consider what processes might 
alter the true shape of the curve. The curve might conceivably be sharpened by 
erosion or removal by slumping. Such effects should be apparent as sudden 
changes in lithology or as time discontinuities. Careful inspection has not re- 
vealed either of these effects associated with the temperature transition in the 
cores of interest. On the other hand, as pointed out by Emiliani (1957), post- 
depositional mixing of the sediment by bottom organisms would tend to spread 
out the change. Such activity is known to occur but is difficult to evaluate 
quantitatively. It may be concluded that the temperature transition is at least 
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Fig. 1. Comparison of surface ocean temperature curves based on foraminifera 
abundance (Ericson and Wollin, 1956a and 1956b; solid curves, warm-cold scale) with 
those based on oxygen isotope measurements (Emiliani, 1955; dotted curves, degrees centi- 


grade scale). 


a 
43) | 
«O° 30° 20° 40° 30° 20° 
cm 
$ 
¢ 
g 
g 
: 100 ; 
¢ 
d 
¢ 
> 
200 ? 9 
~9 
4 J 
° 
300 
» 
400 


Fig. 2. Foraminifera temperature curves and radiocarbon ages for deep sea cores. 
The climate curves were constructed by Ericson (Ericson and Wollin, 1956a and b; Ewing, 
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et al., 1958). The radiocarbon ages on cores A 179-4, A 172-6 and A 180-73 were obtained 
by Rubin and Suess (1955 and 1956). The dates on the other cores were obtained at the 
Lamont Geological Observatory (Broecker, et al., 1956; Broecker and Kulp, 1957). 
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as sharp as indicated by the temperature curves.’ 

Since the carbonate material in the cores of interest consists mainly of 
foraminiferal tests and coccoliths. both of which grow in the surface water, 
C'* age determinations can be made by comparing the C'*/C™ ratio in the core 
material with that in dissolved carbonate in present day surface ocean water. 
Phe difference is attributed to radioactive decay, There are two major poten- 
tial sources of error in such determinations: 1) post-depositional exchange, 2) 
presence of reworked carbonate material, The fact that the oxygen isotope 
record is preserved and that C'' ages of >30,000 years have been observed on 
core materials suggest that the amount of exchange is small and should be 
negligible in the age range of interest. Suess (1956) pointed out that compari- 
son of C'* ages on the coarse and fine fraction should give an indication as to 
the presence of reworked carbonate. In one case he found that the fine fraction 
showed an appreciably greater age. Two such comparisons by Broecker and 
Kulp (1957) showed a greater age for the fine fraction in one case but no 
significant difference in the other. For the age range of interest reworked car- 
bonate appears to be the only important source of systematic error. Thus it 
may be concluded that the ages obtained are either equal to or greater than the 
true age of the core material. 

The micropaleontologic temperature curves obtained by Ericson and 
Wollin (1956) and the radiocarbon dates obtained by Rubin and Suess (1955, 
1956), Broecker, et al. (1956). Broecker and Kulp (1957) are shown in 
figure 2, With the exception of cores A 179-8 and A 185-35, which contain 
numerous turbidity current deposits and core A 180-48 which shows an ab- 
normally high sedimentation rate, the cores may be taken to represent normal 
deep-sea sedimentation. A more complete discussion of the lithology of these 
cores has been published elsewhere (Ericson, et al., 1956; Ericson and Wollin, 
1956a. and 1956b: Ewing, et al., 1958). 

Table 1 summarizes the data. The best estimate of the age of the mid- 
point of the temperature transition in each core is given in column 7, Columns 


5 and 6 give respectively the estimates for the maximum age for the beginning 


of the temperature transition and the minimum age for the end. The difference 
between these two ages provides an upper limit on the duration of the transi- 


tion interval (column 8). 

* Emiliani (1957) concludes that the temperature changes shown in several of the cores 
used in this study (R 10-10, A 180-48, A 179-15) are anomalously sharp and that the more 
gradual changes suggested by some of the cores (A 172-6, A 180-73, \ 180-74) give a bet- 
ter picture of the actual shape of the temperature transition, Since the latter are slow 
deposition cores mixing by bottom borers would be more effective in spreading out the 
time record. Unless erosion removed a portion of the rapid deposition cores at precisely 
the time when the temperature was changing it is dificult to postulate a mechanism for 
sharpening the curves. Neither the lithologie or age data gives any indication of loss of the 
transition period sediment by erosion. For these reasons more weight is given to the cores 
showing sharp transitions. 


Emiliani (1957) has suggested that mixing by bottom borers can make the age at a 
given depth in a core too old, As evidence he uses the finite ages of the tops of several 
cores, However, if the mixing process has proceeded at a more or less constant rate, the 
amount of younger material mixed downward should roughly compensate tor the amount 
of older material mixed upward in samples collected from below the depth of mixing, In 
this case ages at the depths of the temperature break should be unaffected by mixing. 
Also, the existence of rather sharp color changes in many of the equatorial cores suggests 
that the mixing produced by bottom borers is more limited than Emiliani suggests. 
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The majority of the midpoint ages lie close to 11,000 years, Those which 
are significantly older probably reflect the influence of reworked carbonate. Al- 
though the midpoint ages close to 11,000 years may also be somewhat too old, 
the excellent agreement between cores of widely different lithology and geo- 
graphic locations suggests that this effect is small (< 1000 years). 

The maximum time for the transition is in all cases less than 6000 years. 
Since the manner in which the limits were determined, the finite spacing of 
the temperature measurements, and the possibility of homogenization by bot- 
tom borers all tend to increase the magnitude of this range, the smaller times 
(~2000 years) obtained for cores R 10-10 and A 180-48 and the very small 
time (~200 years) for V 12-97 are probably significant.° No minimum esti- 
mate can be made from the data in hand. 

The conclusion is drawn that a large portion of the Atlantic Ocean under- 
went a temperature increase of 6° to 10°C during a period of less than 2000 
years and that the midpoint of this change was within 300 years of 11,000 
years ago. 

Sedimentation rate.—A study of the variation in sedimentation rate with 
time in core A 180-74 from the mid-equatorial Atlantic by Broecker, et al. 
(1958) indicates that both the carbonate and the clay fraction in this core un- 
derwent a significant decrease in sedimentation rate close to 11,000 years ago. 
Cumulative plots of weight of clay and of weight of carbonate as a function of 
time before present as determined by radiocarbon dates (Broecker, et al., 1956; 
Broecker and Kulp, 1957) is shown in figure 3. From the change in slope 
close to 11,000 years ago it can be shown that the clay rate decreased by a 
factor of 3.7 and the carbonate rate by a factor of 2.1. Although this represents 
only one locality and does not permit an extrapolation to other parts of the 
Atlantic, the data supports the hypothesis that an abrupt change in climate 
occurred close to 11,000 years ago. 

Ventilation of restricted basins.—The depletion of oxygen from the bottom 
water in restricted marine basins is often recorded in the sediments by the 
preservation of organic materials, the extermination of the benthic fauna, and 
the presence of H.S gas. Such a situation exists in the upper 5 to 10 meters of 
sediments in the Cariaco Trench off Venezuela, Although the stagnation of the 
deep water in the trench appears to have been continuous during the period of 
deposition of this sediment, the presence of a layer of oxidized clay beneath 
the organic rich ooze demonstrates aeration at some time in the past. A radio- 
carbon date on organic material from the base of the anaerobic layer indicates 
that the abrupt stagnation of the trench occurred within 300 years of 11,000 
years ago (Heezen, et al., 1958 and 1959). 

EVIDENCE FROM FRESH-WATER DEPOSITS 


Deposits formed by fresh-water drainage systems also record the abrupt 
climatic change seen in deep-sea deposits. Two such systems will be considered: 


* Emiliani (1957) concludes that the transition occurred gradually over a period of about 
8000 years. He considers cores R 10-10 and A 179-15 which show rapid transitions un- 
reliable because of the large changes in sedimentation rate suggested by the ages. Ericson 
(personal communication) who did the lithologic studies of these cores, considers them 
both reliable in the depth range of interest. Core A 179-15 shows no evidence of abnormal 
sedimentation above a depth of 300 cm. Even if abnormalities were present it is difficult 
to construct a sequence of events which would yield the observed temperature curves and 
radiocarbon dates if a gradual temperature change over a period of 8000 years did occur. 
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Fig. 3. Cumulative weight of CaCO; (upper diagram) and of clay (lower diagram) 
as a function of time in the past (Broecker, et al., 1958). The points represent radiocarbon 
ages with their respective errors. The slope of the lines connecting points gives the 
sedimentation rate at a given time, For both components a sharp change in rate occurs 
close to 11.000 Vears ago, 


1) the system of lakes in the Great Basin of the western United States and 2) 
the main midcontinent drainage system consisting of the Great Lakes, the 
Mississippi River and the St. Lawrence River. 
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Pluvial lakes.—The Great Basin lakes are of particular interest since they 
are direct indicators of climate, Since these lakes have no outlet, their size is 
controlled primarily by a combination of rainfall and evaporation rates, As 
shown by Broecker and Orr (1958) the lakes respond rapidly to climatic 
change so that their level at a given time provides an excellent index of climate. 
Radiocarbon dates for Searles Lake by Libby (1955) (see Flint and Gale, 
1958) and for lakes Lahontan and Bonneville by Broecker and Orr (1958) 
provide reasonably detailed data on the level of these lakes as a function of 
time. 

Since these curves are based primarily on C** dates obtained on material 
formed in fresh water lakes, the problem of estimating the initial C’* concen- 
tration in such materials must be considered. As shown by Deevey, et al. 
(1954) in some cases the concentration is significantly lower than in terrestrial 
plants or surface ocean water. Broecker and Orr (1958) used the present 
C'/C* ratio in the dissolved carbonate of Pyramid Lake as the control value. 
Broecker and Walton (1959) on the basis of numerous C'* measurements on 
currently forming materials from Great Basin fresh-water systems concluded 
that the C'*/C** ratio in these lakes was within 5 percent of the present Pyra- 
mid Lake value during high water periods, This corresponds to a 400-year 
age uncertainty. 

A second possibility of systematic error, especially for carbonate samples, 
is the incorporation of C'* by the sample after deposition, This could occur 
through transport and redeposition by rain or ground water or by exchange 
with the atmosphere, Experiments by Broecker and Orr (1958) suggest that 
contamination by atmospheric exchange is negligible in nearly all cases, Only 
the dry climate of the Great Basin and the internal consistency of the data can 
be used as evidence against the importance of redeposition, The ages given by 
Broecker and Orr on carbonate materials are probably good to 400 years but 
in some cases may represent only minimum estimates. 

The available radiocarbon data, summarized in tables 2 and 3, strongly 
suggest an abrupt decrease in level from near maximum to near desiccation 
close to 11,000 years ago. The evidence may be summarized as follows, Except 
for the samples from the Winnemucca Cave area, the tufa samples from the 
highest level of Lake Lahontan yielded ages between 11,000 and 12,000 years 
B.P. averaging 11,600 years B.p. Samples from the Winnemucca Cave area 
suggest a second maximum close to 10,000 years ago. Although the Bonneville 
samples show evidence of high lake levels at both of these times no distinct 
separation into two groups exists. The Searles Lake data indicate that the 
now-dry lake contained water prior to 11,000 years ago. Data from each of the 
three lakes clearly demonstrates the return of arid conditions close to 11,000 
years ago. Searles Lake became completely desiccated. Lake Lahontan fell be- 
low the level of Fishbone Cave and Leonard Rock Shelter and Lake Bonneville 
fell below the level of Danger Cave and deposited a sequence of carbonates on 
the salt flats. 

The evidence for a rapid rise following the low levels of these lakes is 
problematic. In both the Lahontan and the Bonneville basins C** data suggest 
a maximum at about 10,000 years ago. No evidence for this maximum appears 
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to be present in the Searles Lake record, however. More work is needed to settle 
this question. 
The events recorded in the deposits of the Great Basin lakes are sum- 
marized in figure 4, A period of continuous relatively high-water levels be- 
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Fig. 4. The height of the major Great Basin lakes as a function of time as interpreted 
from radiocarbon ages. 


ginning about 25.000 years was culminated by a sharp drop in lake level close 
to 11,000 years ago. This drop may have been followed by a brief high-water 
period ending about 9500 years ago. Since that time the lakes have had levels 
similar to those observed at present. 

Midcontinent drainage system.—The midcontinent drainage system re- 
cords climate in a somewhat more indirect manner, The sediments in the al- 
luvial valley of the Mississippi River, as well as those on the Sigsbee Abyssa 
Plain of the Gulf of Mexico. record a marked change in the operation of the 
drainage system. Ewing. et al. (1958) have shown on the basis of sediment 
core studies that great thicknesses of gray silty clay were deposited on the 
abyssal plain during late Wisconsin time. They showed that this silty clay 
emanated from the mouth of the Mississippi River and was transported by 
turbidity currents. Radiocarbon dates (Broecker and Kulp, 1957) on globi- 
gerina ooze from immediately above the silty clay sequence from one of the 
cores (A 185-35) indicates that the supply of silty clay was cut off abruptly 
close to 11,000 years ago. 

The evidence used by Ewing, et al. (1958) to establish the origin of the 
silty clay deposits is summarized in an idealized diagram in figure 5. Cores 
(see, for example, V 3-127, fig. 2) on knolls rising above the abyssal plain 
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Fig. 5. Idealized diagram showing relationships of sediments on the floor of the 
Gulf of Mexico to those in the alluvial valley of the Mississippi River, The Gulf sediments 
are described by Ewing, et al. (1958) and the alluvial valley sediments by Fisk and 
McFarlan (1955). 


are composed entirely of foraminiferal ooze, show normal climatic curves and 
have, as shown by radiocarbon dates, sedimentation rates similar to those in 
other normal deep sea oozes, Since the cores taken on the abyssal plain sur- 
rounding the knolls all show a thick silty clay sequence, Ewing, et al. conclude 
that the sediment must have been carried in along the bottom rather than 
having fallen from the surface. 

Further evidence for an abrupt change in the operation of the Mississippi 
River system is found in the sediments on the alluvial plain in Louisiana, Fisk 
(1955) has classified the sediments filling the channel cut by the river during 
the last period of sea level lowering into two categories; the substratum con- 
sisting of sand and gravel deposits and the top stratum consisting of silt and 
clay. A number of radiocarbon dates (Broecker, et al., 1956; Broecker and 
Kulp. 1957; Brannon, et al., 1957) on samples taken near this transition zone 
indicate that it has an age close to 11,000 years. The results are summarized 
in table 4, It seems reasonable to believe that whereas’ much of the clay and 
silt load of the Mississippi River was swept beyond its channel and beyond the 
continental shelf into the Gulf of Mexico during the late Wisconsin, the char- 
acteristics of the river changed close to 11,000 years ago in such a manner that 
these clays and silts were deposited in the present alluvial valley and in the 
modern delta. In the absence of other evidence it is difficult to say whether this 
change was brought about by a change in sea level, discharge, or quantity and 
character of the detrital load. 


The level of Lake Michigan provides an indirect index of climate, High 
levels are thought to correspond to times when glaciers covered the Straits of 
Mackinac causing the lake to drain to the south, On the other hand, when the 
ice receded to a position north of the glacially depressed straits drainage to 
the east (perhaps via the St. Lawrence) allowed the lake level to fall, Thus 
levels of Lake Michigan at or below its present level indicate periods of rela- 
tively warm climate. 
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Evidence for low water intervals in Lake Michigan close to 11,000 years 
ago is provided by both the Two Creeks Forest dated at 11,300 years B.P. (see 
table 5) and the peat bed exposed at South Haven dated at 11,200+600 
(Crane, 1956). 10.860+350 (Rubin and Suess, 1955), and 10,790+200 
(Preston, et al.. 1955). The use of this data to support a warming in climate 
close to 11.000 years ago is complicated by the Valders readvance, If the age 
averages of 11.300 for Two Creeks and 10,900 for South Haven are accepted 
at face value and if the Two Creeks Forest was overrun by Valders ice shortly 
after growth, then two low levels for Lake Michigan, one preceding and one 
following the Valders. are required. The consequent high advance and retreat 
rates of the Valders ice can be avoided by assuming that the Two Creeks 
Forest and South Haven peat are both pre-Valders in age. This would require 
that either or both of the radiocarbon averages are slightly in error or that a 
significant lag occurred between the growth and destruction of the Two Creeks 
Forest. Again, as in the case of the Great Basin, more work is needed to re- 
solve the multiple events which occurred in the time interval between 10,000 
and 12.000 years ago. 

Radiocarbon ages on marine shell from the St. Lawrence Valley provide 
additional evidence for a rapid ice retreat. In order for these glacially de 
pressed lowlands to be invaded by marine waters the ice must have retreated 
considerably from its maximum extent (see Flint, 1956), The results are as 
follows: Hull, Ontario, Y 215, 10,630 + 330, Ottawa, Ontario, Y 2106. 
10.850 + 330, and Montreal, Quebec, 11,370 + 360 years (Preston. et al.. 
1955). Again the weighted average lies slightly below 11.000 years. 

Leverett and Taylor (1915). however, correlated the Champlain Sea with 
the Lake Algonquin stage of Lake Michigan and more recently Flint (1950) 
correlated it with the Lake Chippewa stage. In either case the Champlain Sea 
invasion would be a much more recent event (4000 to 5000 years B.p.). A 
radiocarbon date of 9430+ 250 on a peat sample from sediments post-dating 
the Champlain Sea (Terasmae, 1959) verifies the antiquity of the shell ruling 
out the explanation given by Flint (1956) that the shell contained considerable 
carbonate derived from ancient limestones. 
EVIDENCE 


FROM POLLEN 


PROFILES 

Another index of climate is the type of plant and animal remains found 
in terrestrial deposits. Of these the pollen profiles of northwestern Europe have 
been most thoroughly studied and are best understood, Also, a large number 
of radiocarbon measurements are available allowing a fairly detailed chronol- 
ogy to be established (see table 6). A portion of this profile is reproduced in 
figure 6. The outstanding feature of all the diagrams is the pronounced warm 
period referred to as Alleréd. It is preceded and followed by colder climates, 
the Older Dryas and Younger Dryas respectively, The available radiocarbon 
data are summarized in table 6. From these average ages for the various pollen 
zones and for the transition can be computed, These averages are shown in 
figure 6. It is clear that the warm Alleréd climates prevailed during the period 
of increase in surface ocean temperature. A problem exists in that there is not 
a unidirectional change in climate from cold to warm following the Older 
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Fig. 6. Sequence of climate periods in the standard European pollen profile (based 
on diagram given by Barendsen, et al., 1957). The average ages for the periods and transi- 
tions are based on the data summarized in table 6, Warm periods are indicated by W and 
cold by C. No radiocarbon data are available for the dotted transitions. 


Dryas for the cold climate is re-established during the period from approxi- 
mately 10,650 to approximately 9800 years ago. The significance of this oscil- 
lation is not clear at present. It is interesting to note that an identical 
oscillation appears in the Great Basin and Great Lakes chronologies. It is en- 
tirely possible that the sharp change in oceanic conditions correlates with the 
post Younger Dryas rather than the pre Younger Dryas warm period as has 
been implied above. 
CONCLUSIONS 

From the evidence listed above it is clear that a major fluctuation in 
climate occurred close to 11,000 years ago, The primary observation that both 
surface ocean temperatures and deep sea sedimentation rates were abruptly 
altered at this time is supplemented by evidence from more local systems. The 
level of the Great Basin lakes fell from the highest terraces to a position close 
to that observed at present. The silt and clay load of the Mississippi River was 
suddenly retained in the alluvial valley and delta, A rapid ice retreat opened 
the northern drainage system of the Great Lakes and terrestrial temperatures 
rose to nearly interglacial levels in northwestern Europe. In each case the 
transition is the most obvious feature of the entire record. 
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Correlation of these systems on a finer scale is difficult in that oscillations 
preceding and following the major change in the case of the European pollen 
profiles, Great Basin’ lakes and Great Lakes region have not been observed in 
the record left in the Mississippi alluvial deposits or deep ocean deposits. More 
work is needed before a satisfactory explanation can be given for these dif- 
ferences. 

Certainly if the abruptness of such a major change in world climate can 
be firmly established it will be of the utmost importance that acceptable 
theories are able to account for it. It is interesting to note that the recent 
“Theory of Ice Ages” proposed by Ewing and Donn (1956 and 1958) was 
stimulated by the observation that the change in climate which occurred at the 
close of the Wisconsin glacial period was extremely abrupt. 
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TABLE 2 
Materials from deposits formed during high water periods 
of the Great Basin lakes 
Height Above 
1890 Lake Sample Age 
Level (feet) Location Material No. (years) 


Pyramid Lake Basin 
Lake Winnemucca tufa L 289.G* 9700+ 200 
Cave Area, Nev. 
Lake Winnemucca tufa L, 356-G* 10,000 + 220 
Cave Area, Nev. 
Lake Winnemucca tufa L 356-H* 9700+ 200 
Cave Area, Nev. 
Lake Winnemucca tufa L 364-AA* 9500 + 200 
Cave Area, Nev. 
Anaho Island tufa L, 289-N* 11,800 + 200 
Pyramid Lake, Nev. 
Anaho Island tufa 289-M * 11,700 + 200 
Pyramid Lake, Nev. 
Anaho Island tufa , 289.L* 11,570+ 25 
Pyramid Lake, Nev. 
Mullen Pass tufa L 289.1* 11,250 +350 
Pyramid Lake, Nev. 
Fishbone Cave tufa L 289-C* 11,700+ 500 
Lake Winnemucca, Nev. 
Anaho Island tufa W 442** 12,050 + 400 
Pyramid Lake, Nev. 
Great Salt Lake Basin 
South end Oquirrh tufa L. 363-D* 11,000 + 600 
Mins., Utah 
South end Oquirrh tufa L 435-G* 11,300 + 250 
Mins., Utah 
South end Oquirrh tufa L 435-C* 10,600 + 300 
Mitns., Utah 
South end Oquirrh tufa W 409** 11,300 + 300 
Mins., Utah 
Pleasant View tufa W 456** 11.650 450 
Salient, Ogden, Utah 
Weber River Delta marl W 382** 12,960 + 3507 
Ogden, Utah 
Weber River Delta peat W 326** 9925 + 3007 
Ogden, Utah (from W 382) 
Weber River Delta peat W 440** 10,260 + 300T 
Ogden, Utah (repeat W 326) 
Searles Lake Basin 
Searles dry lake, organic material C 894*** 10,494. + 560 
California from top of lake 
muds 
Searles dry lake, carbonate from Y 574-a**** 11,810+ 140 
California top of lake muds 
Searles dry lake, organic material Y 574b**** 10,700+ 130 
California from top of lake 
muds 
Broecker and Orr (1958) *** Libby (1955) 
Rubin and Alexander (1958) **** Flint and Gale (1958) 


A new set of samples collected from this locality by one of the authors (Wallace S. 
Broecker) yielded ages of 9500+200 for the organic fraction and 9700+200 for 
the carbonate fraction. The high age obtained by the Washington laboratory for the 
original marl has been checked and confirmed (Rubin, personal communication) 
strongly suggesting that this sample contained reworked detrital carbonate, Since 
these peats are overlain by lake sediments, these dates provide excellent evidence 
for a post 11,000 year p.p, high lake stand. 
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TABLE 3 


Materials from deposits formed during low water periods 


of the Great Basin lakes 


Sample Age 


Mater a No. (vears)} 


Pyramid Lake Basin 
Fishbone Cave twigs overlying lake 11,200 + 
Lake Winnemucca, Nev. sediments 
Leonard Rock Shelter guano overlying , 599* 11,200 + 
Humboldt Lake, Nev. lake sediments 


Needles Area outermost tufa 8500 + 
Pyramid Lake, Nev layer on dome 


Great Salt Lake Basin 
Danger Cave, Bonneville sheep dung over- M 118**** 11,000 + 
Salt Flats, Utah lying beach gravels 
Danger Cave, Bonneville wood overlying C610** 11,150 +57 
Salt Flats, l tah lake de posits 
Danger Cave, Bonneville twigs associated M 119**** 10,400 + 
Salt Flats, Utah with M 118 
Hooper, Utah plant stems in W 386*** 9730+ 3 


growth position 


Bonneville Salt Flats carbonate from one L 485* 11,150 + 
Knolls, Utah foot beneath surface 


™ arles Lake Basin 


Searles ake organic material RC-36***** 10,270+45 
California from mud within 
upper salt 


Searles dry lake organic material RC-50***** 11,400+600 


California from mud within 


Searles dry lake upper salt 
California carbonate from RC-50c**929 9900 + 500 
RC-50 


Searles dry lake trona from base W-248* ** 8550 + 


California of upper salt 


* Broecker and Orr (1958) 
Libby (1955) 
Rubin and Alexander (1958) 
**** Crane (1956) 


***** Flint and Gale (1958) 
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Evidence for an Abrupt Change in Climate Close to 11,000 Years Ago 
TABLE 4 


Radiocarbon dates defining the transition from sand to clay deposition 
in the Mississippi Alluvial Valley 


Depth to 
Sample Transition 
Depth Zone Sample Age 
Location (feet) (feet) No. (years) 


La Fourche Parish, La, ~ 5 HOR-73*** 8800 ++ 180 
Duck Lake Area, La. : L-175B* 9750 + 550 
Dulac Area, La. ‘ : L-291X** 10,700 + 150 
Grand Isle Area, La. 21! 24 L-291G** 11,050 -+ 300 
Grand Isle Area, La, 24 2: L-291H** 10,530 + 350 
Bayou Pigeon Area, La. 220 HOR-126* ** 13,650 + 300 


* Broecker, et al., 1956 
** Broecker and Kulp, 1957 


*** Brannon, et al., 1957 


TABLE 5 


Radiocarbon dates on Two Creeks forest samples 


Laboratory No. of Samples Average Age (years) 


Chicago 11,400+ 350 (Libby, 1955) 
Yale 11,130+350 (Preston et al., 1955) 
Washington ; 11,375+150 (Suess, 1954) 
Michigan 10,550+450 (Crane, 1956) 


Weighted Average  11,300+130 


TABLE 6 
Radiocarbon dates for northeastern Europe pollen profiles 
Sample No. Age (years) 


Younger Dryas 


GRO 419 10,345 + 275 
K 111 10,300 + 350 
Y 139-3 11,350 + 150* 


Younger Dryas/Alleréd Boundary 


10,890 + 240 
10,995 + 250 
10,260 + 200 
10,510 + 180 
10,500 + 400 
10,560 + 200 
H 105-87 11,500 + 300* 
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TasLe 6 (Continued) 
Allerod 


103 11,060 + 480 
113 10,930 + 380 
104 10,990 + 240 
110 10,770 + 360 
105 11,800 + 410 
106 11,880 + 340 
112 11,700 + 360 
107 11,160 + 320 
337 11,044 + 500 
909 11,025 + 120 
907-937 10,650 + 120 
9?0 10,660 + 90 
647 11,230 + 400 
925 11,065 + 120 
948-933 11,700 + 90 
11,470+90 
11,560 + 100 
11,200 + 320 
11,560 + 260 
11,900 + 500 
11,800 + 300 
11,550 + 280 
11,930 + 200 
11,450 + 180 
12,500 + 180* 
10.880 + 160 
11,220 + 350 


Allerod/Older Dryas 
Older Dryas 
Older Drvas/Bolling 
GRO 926 825 + 120 
H 77-54 ey + 260 
Bolling 
GRO 927 2,355+170 
GRO 1104 12,300 + 100 
Bolling/Oldest Dryas 
GRO 928 + 100 
GRO 935 + 130 
H 88-74 280 
H 106-89 320 
H—Hiedelberg (Munnich, 1957) 
GRO—Groningen (DeVries, et al., 1958) 
K—Copenhagen (Anderson, et al., 1953) 
W—Washington (Rubin and Suess, 1955 and 1956) 
Y—Yale (Barendsen, et al., 1957) 


Not included in average in figure 6. 
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